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ABSTRACT  OF  THE  THESIS 


The  Effects  of  a  Geomagnetic 
Storm  on  Thermospheric  Circulation 

by 

Douglas  George  Brinkman 
Master  of  Science  in  Atmospheric  Science 
University  of  California,  Los  Angeles,  1987 
Professor  S.  V.  Venkatesvaran,  Chair 

The  motions  of  the  thermosphere  and  its  interactions  with  the  iono¬ 
sphere  during  a  geomagnetic  storm  are  of  current  interest  to  space 
scientists.  A  two-dimensional  model  was  used  to  simulate  the  thermo¬ 
spheric  response  to  the  Impulsive  high-latitude  heating  associated  with 
a  geomagnetic  storm.  The  storm-induced  motions  can  be  characterized  by 
an  initial  period  of  transient  waves  followed  by  the  development  of  a 
mean  circulation.  These  motions  generate  an  electrical  current  system 
that  is  on  the  same  order  of  magnitude  as,  and  in  the  opposite  sense 
to  the  normal  Sq  current  system.  Model  simulated  winds  and  electrical 
currents  were  then  compared  to  observations. 


XI 


The  motions  of  the  thermosphere  and  Its  Interaction  with  the  iono¬ 


sphere  during  a  geomagnetic  storm  are  of  current  interest  to  space  sci¬ 
entists.  The  thermosphere  is  that  part  of  the  earth's  atmosphere  that 
extends  from  80  km  to  450  km.  It  is  composed  mainlv  of  neutral  partic¬ 
les  with  temperatures  ranging  from  180K  at  the  mesopause  to  well  over 
1000K  at  its  upper  reaches.  An  Integral  part  of  the  thermosphere  is 
the  ionosphere  which  is  composed  of  several  1 avers  of  ionized  partic¬ 
les  that  reside  within  the  thermosphere.  The  response  of  the  thermo¬ 
sphere  to  various  Inputs  and  the  subsequent  interaction  of  neutrals 
and  ionized  particles  can  yield  interesting  results. 

Of  particular  interest  is  how  these  interactions  Influence  long 
range  radio  communications  and  satellite  drag.  The  feasiblitv  and 
quality  of  radio  communications  is  dependent  on  the  state  of  the  iono¬ 
sphere  and  hence  the  motions  of  the  thermosphere.  An  increase  in  drag 
can  greatly  affect  the  orbit  of  low-to-medlum  altitude  satellites.  A 
vivid  example  of  this  was  the  orbit  decay  and  reentry  of  Skvlab. 

The  purpose  of  this  studv  was  to  aquaint  myself  with  the  two  di¬ 
mensional  thermospheric  model  of  Richmond  and  Matsushita  <197S>.  to 
simulate  the  hvdrodynamical  and  electrodynamical  responses  to  a  gene¬ 
ric  geomagnetic  storm  ,  and  to  simulate  the  responses  for  two  actual 
geomagnetic  storms  to  see  how  well  the  model  reproduced  the  time  de¬ 
pendent  storm-time  features  of  the  thermosphere.  Chapter  two  describes 
the  quiet  and  disturbed  thermosphere  including  basic  structure,  energy 
sources,  observation  methods,  and  observations  of  the  quiet-time  and 


storm-time  neutral  wind  circulation.  Chapter  three  describes  the  two- 
dimensional  time  dependent  thermospheric  model  of  Richmond  and 
Matsushita  <1975)  used  in  this  studv.  Basic  characteristics  of  the 
model  and  mode  Inputs  are  discussed.  Chapter  four  discusses  the 
model  outputs  as  they  apply  to  a  generic  geomagnetic  storm.  Chapter 
five  describes  some  of  the  electrodvnamic  aspects  of  the  thermosphere 
and  how  it  is  affected  by  the  storm-time  circulation.  Chapter  six  com¬ 
pares  two  geomagnetic  storms  analyzed  bv  Mazaudier  et  al .  (1985)  with 
model  outputs.  Finallv,  chapter  seven  provides  conclusions  and  dis¬ 
cusses  plans  for  future  studies  of  thermospheric  processes. 


II  BASIC  PROPERTIES  OF  THE  QUIET  AND  DISTURBED  THERMOSPHERE 
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The  structure  of  the  thermosphere  is  the  result  of  a  complex  Inter¬ 
action  between  solar  energy  Inputs,  composition,  and  circulation.  It 
is  difficult  to  describe  anv  one  part  of  this  system  without  referring 
to  the  effects  of  the  other  parts.  Several  excellent  reviews  of  the 
thermosphere-ionosphere  system,  have  been  written  bv  Hanson  and  Carlson 
<1977).  Killeen  <1987),  Richmond  <1983),  Rishbeth  <1972),  and  Roble 
(1977) . 

This  section  brieflv  summarizes  that  material  and  is  not  meant  to 
be  an  all  inclusive  discussion  of  the  thermosphere-ionosphere  system. 
The  effects  of  solar  radiation  Inputs  on  composition  and  temperature 
are  described.  Energv  inputs  from  the  solar  wind  are  discussed.  This 
is  followed  bv  a  description  of  the  resulting  circulation  from  both 
these  energy  Inputs.  To  round  out  the  discussion,  observational  tech¬ 
niques  are  described  and  a  summary  of  neutral  wind  observations  for 
both  quiet  and  storm-times  are  presented. 

Energy  from  the  sun  can  be  divided  into  two  distinct  forms,  that 
which  is  associated  with  electromagnetic  radiation  and  that  which  is 
associated  with  the  solar  wind.  Most  of  the  sun’s  electromagnetic 
radiation  is  in  the  visible  part  of  the  spectrum  and  passes  through  the 
thermosphere  to  the  ground  and  lower  atmosphere  where  it  is  absorbed 
or  reflected  back  to  space.  Solar  electromagnetic  radiation  with  wave¬ 
lengths  less  then  200nm  is  absorbed  by  constituents  in  the  thermosphere 


such  as  0,  0^,  ,  and  NO. 
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Composition 


Below  80  Km.  the  composition  of  the  the  atmosphere  Is  7855  2135 
O^.  and  1 1  of  the  other  constituents  such  as  0^,  CO^,  and  H^O.  Neutral 
gas  composition  above  80  km  Is  dependent  on  solar  extreme  ultraviolet 
<EUV>  radiation  and  ultraviolet  <UV>  radiation  with  wavelengths  less 
than  200nm  which  photodlssoclates  molecules  Into  their  atomic  form. 
The  reaction: 

C>2  ♦  hv  =  0  +  0 

Is  verv  efficient  and  is  the  dominate  reaction  at  these  levels.  Recom¬ 
bination  of  oxvgen  atoms  back  to  an  oxvgen  molecule  requires  a  third 
bodv  such  that: 

o  +  o+  m  =  o2  +  m. 

This  reaction  Is  very  slow  above  90  km.  Oxygen  atoms  must  therefore, 
be  transported  downward  In  order  to  recombine.  This  leads  to  oxygen 
atoms  being  the  major  constituent  of  the  thermosphere  at  150  km. 

Contributing  to  this  maximum  of  atomic  oxygen  is  the  turbopause, 
which  is  at  roughly  110  km.  Above  this  level  molecular  diffusion  be¬ 
comes  important.  As  a  result  constituents  become  segregated  In  alti¬ 
tude  according  to  molecular  weight.  Heavier  species  maximize  at  lower 
themospherlc  altitudes  while  lighter  species  such  as  0  become  dominant 
at  higher  altitudes.  Diffusive  equilibrium  Is  dependent  on  tempera- 


ture,  and  strong  diurnal  variation  of  temperature  causes  the  composit¬ 
ional  structure  of  the  thermosphere  to  exhibit  a  corresponding  strong 
diurnal  variation. 

Temperature 

Between  ISO  and  300  km.  neutral  heating  through  the  absorption  of 
solar  EUV  by  oxygen  atoms  dominates.  Above  300  km,  EUV  absorption  is 
stil  1  important  but  the  effects  of  magnetospherlc  processes  also 
matter.  Below  ISO  km,  absorption  of  solar  UV  In  the  Schumann-Runge 
contlnumun  <130-17Snm>  dominates,  and  below  100  km  0^  absorption  in  the 
Schumann-Runge  bands  (17S-200nm>  Is  important.  Below  80  km,  0„  becomes 

J 

the  dominant  absorber. 

As  described  above,  there  are  a  wide  variety  of  absorbers  In  the 
thermosphere.  However,  above  90  km  there  are  few  effective  radiators 
of  thermal  energy.  During  solar  minimum.  The  global  mean  vertical 
temperature  structure  Is  determined  by  a  balance  between  local  heating, 
and  molecular  and  eddy  thermal  conduction  downward  to  the  mesopause  at 
80  km  where  CO^  and  Og  can  radiate  this  energy  back  to  space.  The  re¬ 
sultant  temperature  ranges  from  180-200K  at  80  km  to  a  temperature 
that  reaches  to  1500K  and  above  at  300  km.  Molecular  conduction  at  the 
upper  levels  Is  so  fast  that  large  vertical  temperature  gradients  can¬ 
not  form  and  the  thermosphere  becomes  isothermal.  The  temperature  at 
this  level  Is  know  as  the  exospheric  temperature.  (The  exosphere  being 
that  part  of  the  atmosphere  where  collisions  are  so  rare,  that  partic¬ 
les  move  In  ballistic  orbits  constrained  only  by  gravity. ) 


The  level  of  solar  activity  also  plays  a  role  In  determining  temp¬ 
erature  structure.  It  Is  associated  with  two  time  scales.  One  time 
scale  la  the  27  dav  solar  rotation  period  and  the  other  Is  the  11  vear 
sunspot  cycle.  Active  regions  such  as  sunspots,  plages,  and  coronal 
condensations,  which  are  visible  on  the  solar  disk,  rotate  with  the 
sun.  They  are  more  prevalent  at  maximums  In  the  11  year  solar  cycle 
and  are  Indicative  of  enhanced  solar  output.  This  enhancement  causes 
the  solar  EUV  radiation  to  double  from  solar  minimum  to  solar  maximum, 
resulting  in  a  TOOK  change  In  the  mean  exospheric  temperature  from  mi¬ 
nimum  to  maximum.  Also  during  solar  max,  radiative  cooling  due  to  HO 
can  play  a  significant  role  in  determining  thermospheric  temperature. 
In  addition,  the  exospheric  temperature  varies  greatly  from  day  to 
day.  Figure  1.  shows  exospheric  temperatures  for  different  levels 
of  geomagnetic  activity. 

Temperature  maximizes  at  the  subsolar  point  (with  a  slight  lag 
due  to  conduction  effects)  and  minimizes  at  the  antisolar  point.  On  a 
seasonal  basis,  the  temperature  maximum  follows  the  subsolar  point 
through  its  yearly  cycle.  Maximum  temperatures  are  observed  at  the 
equinoxes  and  minima  at  solstices.  The  diurnal  temperature  variation 
Is  on  the  order  of  20-302  of  the  mean  exospheric  temperature.  This 
large  diurnal  variation  demonstrates  the  high  degree  of  control  that 
solar  electromagnetic  radiation  exercises  In  determining  thermospheric 
temperature . 


6 


Ionosphere 


The  earth  is  surrounded  by  a  layer  of  ionized  gas  called  the  iono¬ 
sphere  which  resides  within  the  thermosphere.  The  ionosphere  can  be 
divided  into  two  distinct  parts:  a  low  latitude  ionosphere  dominated 
bv  solar  radiation  (which  is  discussed  here)  and  a  polar  ionosphere 
dominated  by  energetic  particle  precipitation  and  magnetospherlc  pro¬ 
cesses  (which  is  discussed  in  the  following  section). 

Several  distinct  layers  of  charged  particles  exist.  The  D  region 

2  4-3 

extends  from  60-90  km  (ion  density  of  10  -10  cm  during  the  day).  The 

5  -3 

E  region,  which  extends  from  90-140  km  (10  cm  ),  is  the  result  of 

the  absorption  of  soft  solar  X-rays.  The  F  region,  which  is  divided 

into  the  FI  layer  from  140-200  km  and  the  F2  layer  from  200  km  to  its 
6  —3 

peak  at  300  km  (10  cm  ),  is  the  result  of  the  absorption  of  solar  EUV 
radiation  from  20-90nm.  Molecular  ions  dominate  in  the  E  and  FI  re¬ 
gions,  while  0+  is  the  dominate  species  in  the  F2  region. 

Chapman  (1931)  first  purposed  that  the  ionization  peaks  were  the 
result  of  the  balance  between  solar  photoionization  and  recombination. 
This  explains  the  E  and  FI  region  peaks  but  not  the  F2  region  peak. 
It  is  downward  ion  diffusion  that  effectively  puts  a  lid  on  the  F2 
laver  and  forms  the  peak. 

There  is  also  considerable  variation  in  ionospheric  density. 

Irregularities  exist  at  all  latitudes  but  predominately  at  high  lati- 

2 

tudes  and  at  the  equator.  Ion  density  esn  vary  on  the  order  of  10 
3  -3 

to  10  cm  within  a  few  kilometers. 


Auroral  Processes 


The  aurora  Is  a  sign  that  large  amounts  of  energy  are  being  depos¬ 
ited  In  the  thermosphere,  especially  at  high  latitudes.  Heating  due 
to  auroral  processes  can  have  a  great  influence  on  the  basic  state  of 
the  thermosphere  and  Its  circulation.  During  quiet  times,  the  high 
latitude  heating  rate  is  about  1/10  of  the  solar  heating  rate  <SX10*^ 
vatts  compared  to  4X10^  watts  at  120  km). 

Energy  from  the  solar  wind  Is  transferred  to  the  neutral  atmosphere 
through  energetic  particle  precipitation,  frictional  dissipation  of 
Ionospheric  currents,  and  Ion  drag  forces,  all  of  which  result  from 
the  Interaction  of  the  solar  wind  with  the  earth's  magnetic  field. 

Particle  precipitation,  as  evidenced  by  the  aurora,  is  the  result 
of  10-0.5  Kev  electrons  traveling  along  field  lines  and  colliding  with 
oxvgen  atoms  In  the  thermosphere.  This  leads  to  the  following  react¬ 
ion: 


O  ♦  e  =  O  ♦  2e  (Ionization) 

The  more  energetic  the  electrons  the  deeper  down  Into  the  atmosphere 
their  energv  is  deposited.  The  resulting  ionization  provides  a  horse¬ 
shoe  shaped  region  of  high  electrical  conductivity  known  as  the  auroral 
oval . 

Solar  wind  Interaction  with  the  earth's  magnetic  field  can  drive 
currents  along  magnetic  field  lines.  Dlscontlnultes  in  the  magneto- 
spheric  convection  electric  field  (^  2*0)  can  lead  to  field-aligned 


potential  differences  that  generate  field-aligned  currents  (  Lyons, 
1980).  The  field-aligned  currents  close  through  the  auroral  oval, 
where  the  frictional  dissipation  of  these  currents  results  In  Joule 
heating . 

Solar  wind  Interaction  with  the  magnetosphere  also  results  In  the 
generation  of  an  electric  field  perpendicular  to  the  magnetic  field. 
High  conductivity  along  field  lines  allows  the  electric  field  to  be 
mapped  down  to  the  Ionosphere.  This  provides  a  momentum  source  that 
sets  the  Ions  and  the  electrons  Into  motion.  Electric  fields  at  high 
latitudes  range  between  10  and  100  mV/m,  which  can  Impart  electron 
drifts  on  the  order  of  200  to  2000  m/s.  The  higher  values  are  seen 
during  intense  geomagnetic  activity. 

The  Ions  collide  with  neutrals  and  drag  the  neutrals  along  their 
drift  paths  <  ion  drag).  This  is  an  important  source  of  momentum  at 
high  latitudes.  The  collisions  between  ions  and  neutrals  also  provide 
heating  and  act  as  mechanical  resistance  to  neutral  wind  motion. 
Hence,  ion  drag  provides  both  a  heat  source  and  a  momentum  source  to 
the  thermosphere,  although  the  effects  tend  to  be  confined  to  the  high 
latitude  region. 

A  geomagnetic  storm  can  greatly  enhance  high  latitude  heating  to  a 
level  equal  to  or  exceeding  heating  from  solar  EUV  and  UV.  By  a  geo¬ 
magnetic  storm,  what  is  meant  is  a  disturbance  in  the  solar  wind  that 
increases  interaction  with  the  earth's  magnetic  field.  It  results  in 
the  development  of  large  electric  fields  and  intense  electrical  cur¬ 
rents  enhancing  high  latitude  heating.  A  geomagnetic  storm  can  last 
from  several  hours  to  several  days  and  show  up  as  a  disturbance  in 
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ground-based  magnetic  measurements.  The  source  of  the  solar  wind  dis¬ 
turbance  could  be  a  solar  flare,  a  coronal  hole,  or  solar  mass  eject¬ 
ion.  By  enhancing  the  auroral  processes  already  discussed,  a  geomag¬ 
netic  storm  can  cause  large  variations  in  compostion,  temperature, 
and  circulation.  These  effects  lag  the  storm  by  several  hours  and  the 
lag  time  is  longer  the  further  one  is  away  from  the  auroral  zone.  The 
amplitude  of  these  changes  depends  on  the  magnitude  of  the  storm. 
Energy  deposited  in  the  auroral  zone  is  transferred  to  the  rest  of  the 
globe  via  dynamic  processes,  and  results  in  a  global  increase  in  exo¬ 
spheric  temperature. 

Circulation 

One  of  the  earliest  ways  of  getting  a  handle  on  thermospheric  cir¬ 
culation  was  through  the  use  or  semi-empirical  models.  The  influence 
of  atmospheric  drag  on  earth  orbiting  satellites  is  a  useful  tool  to 
study  the  thermosphere.  Drag  is  related  to  atmospheric  density. 
Jacchia  (1365)  used  orbital  changes  of  satellites  due  to  drag  effects 
to  deduce  a  global  distribution  of  neutral  gas  density  above  120  km. 
This  density  distribution  was  then  used  to  determine  an  exospheric  tem¬ 
perature  pattern  which  can  be  used  as  an  input  to  a  semi-empirical  mo¬ 
del  such  as  developed  by  Dickinson  et  al .  <1975),  Hedln  et  al.  (1977), 
and  Hedin  (1983).  (Hedin  (1983)  is  an  updated  version  of  Hedin  et  al . ) 
These  models  use  the  temperature  pattern  and  composition  to  derive  a 
pressure  field  that  in  turn  drives  a  wind  field. 

During  quiet  times,  when  solar  EUV  and  UV  are  the  dominant  heat 


sources,  winds  blow  from  the  subsolar  point  to  the  antisolar  point. 
This  circulation  exhibits  generally  upward  motion  on  the  dayslde  and 
downward  motion  on  the  nlghtslde,  with  meridional  winds  blowing  from 
the  equator  toward  the  pole,  and  zonal  winds  blowing  eastward  in  the 
afternoon  and  westward  in  the  morning.  (See  figure  2.,  which  is  from 
Roble,  1977;  but  was  originally  from  Dickinson  and  Geisler,  1968.)  On 
a  seasonal  basis  the  flow  is  from  the  equstor  to  the  poles  at  equinox 
and  from  the  summer  hemisphere  to  the  winter  hemisphere  at  solstice. 

There  can  be  large  departures  from  the  mean  flow  at  high  latitudes, 
where  auroral  heating  is  large  compared  to  solar  heating.  Adding  a 
heat  source  at  high  latitudes  sets  up  a  reverse  circulation  that  coun¬ 
ters  the  quiet-time  winds.  Additionally,  variation  in  the  heating 
source  during  a  geomagnetic  storm  can  generate  gravity  waves.  At  equi¬ 
nox,  the  circulation  due  to  high  latitude  heating  drives  an  average 
flow  from  high  latitudes  to  low  latitudes  in  both  hemispheres  as  shown 
in  figure  3.,  which  is  a  schematic  diagram  from  Roble  (1977),  based  on 
computer  simulations  from  Dickinson  et  al.  (1975).  Its  magnitude  de¬ 
pends  on  the  level  of  geomagnetic  activity.  At  solstice,  Increased 
geomagnetic  activity  enhances  the  flow  in  the  summer  hemisphere  and  re¬ 
verses  the  flow  in  the  winter  hemisphere  as  shown  in  figure  4. (same 
source  as  figure  3.).  On  a  day-to-day  basis,  however,  the  flow  is 
highly  variable. 

Superimposed  on  these  motions  are  the  effects  of  upward  propagat¬ 
ing  planetary  waves,  gravity  waves,  and  tides.  Tidal  oscillations  in 
particular  can  play  a  role  in  lower  thermospheric  motions.  These  os¬ 
cillations  are  excited  by  the  strong  diurnal  variation  of  EUV  and  UV 
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LOCAL  TIME 

WIND  VECTORS  AT  THE  300  KM  LEVEL 


Figure  2.  The  global  distribution  of  neutral  gas  temperatures  at  300 
Kn  based  on  an  empirical  model.  The  arrows  show  the  wind  pattern  de¬ 
termined  at  300  Km.  The  longest  arrow  represents  a  wind  speed  of  22S 
m/s  (from  Dickinson  and  Gelsler,  1968;  and  Roble,  1977). 
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absorption  In  the  thermosphere  and  the  absorption  of  UV  by  ozone  and 
H^O  in  the  mesosphere  and  stratosphere.  Some  tidal  modes  remain  trap¬ 
ped  and  don't  propagate  far  from  their  generation  place.  Other  modes 
propagate  upwards  and  impart  their  energy  to  the  lower  thermosphere. 
The  diurnal  and  semi-diurnal  motions  Induced  by  these  tides  can  gener¬ 
ate  electrical  currents  in  the  ionosphere. 

Other  forces  that  affect  thermospheric  circulation  Include  viscous 
forces  and  ion  drag  forces.  Kinematic  molecular  viscosity  increases 
exponentially  with  height,  and  above  300  km  motions  are  characterized 
by  bulk  motion. 

Ion  drag  also  affects  air  motions.  On  the  dayslde,  ion  density 
maximizes  and  results  in  collisions  decreasing  the  wind  speed.  It  also 
causes  cross-isobarlc  flow  from  high  to  low  pressure  areas.  At  night 
when  ion  density  is  lower  the  effects  of  ion  drag  are  lower  and  wind 
speeds  are  higher.  Quiet-time  circulation  is  hence  a  balance  between 
pressure  gradient,  corlolls,  and  ion  drag  forces. 

Observation  Techniques 

Observations  of  the  thermosphere  are  made  in  several  ways.  The 
earliest  measurements  of  neutral  winds  were  made  by  radar  observations 
of  meteor  trails  at  80-100  km.  This  showed  the  influence  of  planetary 
waves,  atmospheric  tides,  and  gravity  waves  in  the  lower  thermosphere 
and  upper  mesosphere.  A  second  method  was  by  observing  the  drifts  of 
chemical  releases  from  rockets.  These  drifts  could  then  be  used  to  de¬ 
duce  mean  zonal  and  meridional  winds.  As  mentioned  earlier,  the 


effects  of  drag  on  low  and  medium  orbiting  satellites  have  been  used 
to  deduce  exospheric  temperature,  and  hence  thermospheric  winds.  More 
recently,  satellites  In  the  Atmospheric  Explorer  series  and  the  DE 
satellites,  carrying  spectrometers,  have  made  measurements  of  Iono¬ 
spheric  and  thermospheric  properties. 

Ground-based  measurements  can  be  made  by  Incoherent  scatter  radars 
(ISR)  and  Fabry-Perot  Interferometers  (FPI).  Incoherent  scatter  ra¬ 
dars  measure  the  energy  backscattered  from  plasma  density  Irregular¬ 
ities  due  to  Thomson  scattering.  The  doppler  shift  associated  with 
this  returned  energy  provides  a  direct  measure  of  Ion  drifts.  The  mag¬ 
nitude  and  shape  of  the  returned  energy  spectra  provides  Ion  density 
and  plasma  temperature.  These  data  are  then  used  to  Infer  meridional 
neutral  winds.  Most  of  the  ISR  facilities  are  monostatic  (one  antenna) 
and  are  restricted  to  measuring  only  one  component  of  Ion  velocity. 
At  Millstone  Hill,  MA,  for  example,  vertical  ion  drifts  are  measured. 
By  subtracting  out  effects  of  amblpolar  diffusion  and  electric  fields 
(which  are  estimated),  the  meridional  component  of  the  neutral  wind  is 
inferred. 

The  French  quadristatlc  Incoherent  scatter  radar  facility  at  St. 
Santin,  France,  and  more  recently  the  EISCAT  ISR  in  Scandinavia  are 
unique  In  that  they  are  multistatic  radars  (multiple  antennas)  and  can 
give  the  total  Ion  velocity  vector  for  the  E  and  F  regions.  Below  110 
Km,  the  neutral  wind  velocity  Is  approximately  equal  to  the  Ion  veloci¬ 
ty.  Between,  110  km  and  180  km,  the  ion  velocity  parallel  to  the  mag¬ 
netic  field  (Vjj)  is  equal  to  the  parallel  component  of  the  neutral 
wind  (V^jj).  Multiplying  Vj|  by  a  geometrical  factor,  which  Is  depend- 


ent  on  station  location,  the  meridional  velocity  is  obtained.  Above 
180  km  the  effects  of  diffusion  must  be  considered  in  calculating  V^. 
The  electric  field  can  also  be  determined.  In  the  F  region: 

=  <2a  x  3>/B2  . 

Thus  by  measuring  (the  ion  velocity  perpendicular  to  the  magnetic 
field  line),  the  electric  field  can  be  determined.  High  conductivity 
along  magnetic  field  lines  allows  the  F  region  electric  field  to  be 
mapped  down  to  the  E  region. 

The  use  of  incoherent  scatter  radars  is  limited  to  daytime  in  the 
E  region,  when  ion  densities  are  high,  but  at  night  F  region  measure¬ 
ments  can  still  be  made.  Also  at  night,  ground-based  measurements  can 
be  made  using  a  Fabry-Perot  interferometer  <FPI).  Neutral  winds  are 
obtained  by  measuring  the  doppler  shift  of  the  630nm  red  line  emissions 
of  O('D)  at  2S0  km  (the  alrglow  level).  By  looking  at  the  broadening 
of  the  doppler  spectrum,  temperature  is  obtained.  FPIs  are  able  to 
make  measurements  in  both  the  meridional  and  zonal  directions.  The 
majority  of  neutral  wind  measurements  have  been  made  by  this  technique 

Observed  Quiet-Time  Neutral  Winds  and  Temperatures 

The  majority  of  quiet-time  neutral  wind  observations  (summarized 
in  Table  1.)  have  been  made  at  mid-latitudes.  Salah  and  Holt  (1974) 
using  the  incoherent  scatter  radar  at  Millstone  Hill,  MA  <42.6N,71.5W> 
measured  equatorward  neutral  winds  between  100  and  ISO  m/s  during  the 


night  at  300  ka  and  poleward  winds  between  25  and  50  m/s  during  the 
day.  Roble  et  al .  (1974),  again  using  the  Millstone  Hill  ISR,  measur¬ 
ed  neutral  winds  at  various  altitudes  between  2S0  and  480  km.  They, 
also,  measured  equatorward  winds  that  averaged  100-1S0  m/s  during  the 
night  and  poleward  winds  between  25  and  50  m/s  during  the  day.  They 
found  that  the  equatorward  winds  maximized  between  21  and  01  local 
time  <LT>  and  then  decreased  In  magnitude  toward  morning.  The  neutral 
winds  were  also  weaker  during  the  winter  than  during  the  summer. 

Hernandez  and  Roble  (1976)  using  a  Fabry-Perot  Interferometer  at 
Fritz  Peak,  CO  MON,  120W)  looked  at  the  seasonal  variation  of  neutral 
winds  during  quiet  times.  They  found  that  equatorward  winds  and  tem¬ 
peratures  maximized  during  summer  and  were  at  a  minimum  during  winter. 
Also,  the  meridional  winds  were  poleward  during  the  day  and  equatorward 
at  night.  Additionally,  the  zonal  winds  changed  In  direction  from 
eastward  to  westward  around  local  midnight 

Hernandez  and  Roble  (1977)  also  looked  at  seasonal  variations  of 
qulet-tlme  neutral  winds.  They  found  that  the  transition  from  pole- 
ward  winds  to  equatorward  winds  came  earlier  In  the  night  the  closer 
It  was  to  summer.  However,  during  winter,  the  zonal  winds  did  not 
change  from  eastward  to  westward.  Winds  measured  during  Sept.  1975 
show  average  equatorward  winds  of  100-200  m/s  (maximizing  at  midnight) 
and  average  zonal  winds  of  2S-S0  m/s.  During  winter,  the  maximum  equa- 
torVard  winds  decreased  in  magnitude  from  100  m/s  In  October  to  20  m/s 
In  December. 

Slpler  et  al.  (1981)  using  a  Fabry-Perot  Interferometer  at  Laurel 
Ridge,  PA  (48. IN, 79. 2W)  measured  night  time  neutral  winds  for  14  Au- 


gust  1980.  They  found  that  the  neutral  temperature  decreased  from 
1S00K  at  21 LT  to  1200K  at  05LT.  The  zonal  winds  changed  from  SO  m/s 
eastward  at  21LT  to  zero  at  OOLT  and  then  became  westward.  The  meri¬ 
dional  winds  went  to  zero  at  21LT  and  then  became  equatorward  reach¬ 
ing  a  maximum  of  75  m/s  at  OOLT.  The  meridional  winds  then  remained 
constant  until  OSLT. 

Hernandez  and  Roble  (1984)  looked  at  the  changes  In  the  thermo¬ 
spheric  F  region  neutral  winds  from  solar  mlmlmum  to  solar  maximum 
over  Fritz  Peak,  CO.  The  measurements  spanned  the  years  from  1973 
until  1979.  Hernandez  and  Roble  found  that  the  observed  wind  pattern 
changed  from  solar  min  to  solar  max.  First  of  all,  during  winter  for 
both  solar  min  and  max,  the  zonal  winds  blew  eastward  throughout  the 
night,  not  reversing  to  westward  until  morning.  During  summer  at 
solar  min,  the  zonal  winds  would  be  westward  throughout  the  night  ave¬ 
raging  50-100  m/s.  While  at  solar  max,  winds  would  be  eastward  during 
the  early  evening,  shifting  to  westward  around  21-23LT  with  a  maximum 
magnitude  of  100  m/s.  The  transition  time  from  eastward  to  westward 
winds  moved  from  early  morning  in  winter  to  early  evening  In  summer. 

During  solar  max,  there  was  a  difference  between  the  magnitude  of 
the  winds  measured  north  of  the  observing  station,  and  winds  measured 
south  of  the  station.  During  summer,  winds  measured  north  of  the 
station  were  100  m/s  larger  than  winds  measured  south  of  the  station. 
During  winter,  this  difference  decreased  to  50  m/s.  The  north-south 
difference  was  not  apparent  during  solar  min.  The  maximum  equatorward 
winds  during  solar  max  ranged  from  60  m/s  (north  of  the  station)  to 
30  m/s  (south  of  the  station)  in  winter,  and  175  m/s  (north)  to  75  m/s 


(south)  In  summer  (maximum  equatorward  winds  In  August).  This  differ¬ 
ence  was  attributed  to  enhanced  auroral  heating  and  ion  drag,  the  in¬ 
fluence  of  which  diminishes  away  from  high  latitudes.  During  solar  min, 
the  maximum  equatorward  winds  ranged  from  7S  m/s  in  winter  to  ISO  m/s 
in  summer. 

Sipler  et  al .  (1982)  also  looked  at  the  variations  in  neutral  winds 
during  solar  min  and  max,  with  measurements  spanning  the  years  1975  to 
1979.  Their  wind  measurements  generally  agreed  with  those  from  Fritz 
Peak.  In  addition,  they  noticed  variations  in  thermospheric  tempera¬ 
ture.  During  solar  max,  the  temperature  maximized  at  approx.  1500K  in 
the  evening  (18LT)  and  decreased  through  the  night  to  1000K  at  06LT. 
During  solar  min,  this  pattern  varied  depending  on  season.  During  win¬ 
ter,  the  temperature  decreased  from  approx.  700K  at  17LT  to  600K  at 
02LT.  During  equinox,  temperatures  decreased  from  900K  at  19LT  to  600K 
at  22LT.  And  during  summer,  temperatures  decreased  from  HOOK  at  20LT 
to  700K  at  23LT. 

The  quiet-time  neutral  wind  pattern  at  low-latitudes  has  been  in¬ 
vestigated  by  Harper  (1973),  Burnside  et  al .  (1981),  and  Burnside  et 
al.  (1983)  using  a  combination  of  Incoherent  scatter  radar  and  Fabry- 
Perot  Interferometer  measurements  at  Areclbo,  Puerto  Rico  (18N.67W). 
The  observations  were  generally  similar  to  those  at  mid-latitudes  with 
zonal  winds,  eastward  in  the  early  evening,  transitioning  to  westward 
in  the  early  morning  hours  during  sumsser.  During  winter,  the  zonal 
winds  remained  eastward  through  the  night,  but  decreased  in  strength 
towards  morning.  (eg.  winds  eastward  at  12S  m/s  @21LT  decreasing  to 
zero  at  04LT. >  The  meridional  winds  were  typically,  50-100  m/s  equa- 


torward  maximizing  at  21LT.  However,  the  meridional  winds  are  compli¬ 
cated  by  tb'  equatorial  midnite  temperature  maximum.  The  maximum  is 
the  result  of  semi-diurnal  tides  and  leads  to  a  density  bulge  at  the 
equator.  This  results  In  convergence  of  the  normal  equatorward  winds 
and  poleward  winds  from  the  bulge  over  Areclbo  at  02-03LT. 

Havs  et  al .  <1979),  using  a  Fabry-Perot  interferometer  at  College, 
Alaska  < 64 . 8N, 147 . 8V)  measured  high-latitude  F  region  neutral  winds 
during  geomagnetic  quiet-times.  They  reported  100  m/s  westward  winds 
in  the  evening  and  100  m/s  eastward  winds  from  00  to  07LT,  along  with 
equatorward  winds  throughout  the  night  with  a  maximum  of  200  m/s. 
They  attributed  their  measurements  to  ion  drag  forces,  which  caused  the 
neutral  particles  to  follow  the  magnetospheric  convection  pattern. 
Killeen  et  al.  <1986)  using  FPI  measurements  from  the  DE2  low  altitude 
polar  orbiting  spacecraft  and  measurements  from  ground-based  FPIs  also 
determined  that  the  high-latitude,  neutral  winds  are  the  result  of 
neutrals  being  dragged  along  by  ions  (ion  drag)  which  follow  the  large 
double  vortex  pattern  driven  by  magnetospheric  convection. 

In  summary,  mid-latitude  quiet-time  meridional  winds  are  poleward 
during  the  day  ,  and  equatorward  at  night.  Zonal  winds  are  eastward 
in  the  afternoon,  and  westward  in  the  morning.  This  pattern  is  gener¬ 
ally  followed  at  low-latitudes,  except  for  the  effects  of  the  equator¬ 
ial  midnite  temperature  maximum,  which  gives  rise  to  a  poleward  wind 
at  02-03LT  and  converges  with  the  normal  equatorward  winds.  At  high- 
latitudes,  the  neutral  winds  follow  the  magnetospheric  convection 


pattern. 


TABLE:  1.  SUMMARY  OF  QUIET-TIME  OBSERVATIONS 


1.  Meridional  winds:  Poleward  (2S-S0  m/s>  during  the  day. 

Equatorward  (100-1S0  m/a)  at  night. 

2.  Meridional  winds  change  fro*  poleward  to  equatorward,  and 
zonal  winds  change  from  eastward  to  westward  at  approx, 
local  midnight. 

3.  Neutral  temperature  decreases  at  night. 

4.  Nighttime  wind  speeds  are  stronger,  and  changes  In  direction 
are  earlier  In  the  night  during  summer  than  during  winter. 

5.  Neutral  winds  are  stronger,  and  show  greater  latitudinal  var¬ 
iation  during  solar  maximum  than  during  solar  minimum. 

6.  Low-latitude  winds  are  slmlllar  to  mid-latitudes,  but  are 
complicated  by  the  equatorial  midnight  temperature  bulge. 

7.  Hlgh-latltude  neutral  winds  follow  the  magnetospherlc 
convection  pattern. 


Storm  Time  Neutral  Winds  and  Temperatures 


During  active  geomagnetic  periods,  energy  input  at  high  latitudes 
can  chanage  the  neutral  vlnd  circulation  and  thermospheric  temperature 
considerably.  The  general  characteristics  of  this  disturbed  circula¬ 
tion,  described  earlier,  are  supported  by  the  following  observations 
(summarized  in  Table  2>.  Hays  and  Roble  (1971)  measured  neutral  winds 
during  two  geomagnetic  storms  using  a  Fabry-Perot  interferometer.  The 
first  storm  from  Oct  31  to  Nov  1,  1968  was  characterized  by  a  stable 
red  arc.  The  measured  winds  were  equatorward  at  250-300  m/s.  The  se¬ 
cond  storm  on  the  night  of  Nay  14-1S,  1969  showed  a  temperature  in¬ 
crease  from  1400K  at  21LT  to  1600K  at  02LT.  The  measured  winds  were 
equatorward  at  350-400  m/s.  In  both  cases,  the  equatorward  winds  were 
significantly  larger  (by  a  factor  of  2)  than  during  quiet-times. 

Hernandez  and  Roble  (1976b),  at  Fritz  Peak,  looked  at  four  geomag¬ 
netic  storms  during  1974  using  a  FPI.  They  measured  equatorward  winds 
averaging  200  to  300  m/s.  For  one  extremely  intense  storm  <Kp=9>,  they 
measured  an  equatorward  velocity  of  640  m/s.  The  zonal  winds  averaged 
100-200  m/s  westward,  and  didn't  exhibit  the  night-time  transition  from 
east  to  westward  winds  so  often  seen  during  quiet  times.  However,  when 
geomagnetic  activity  decreased,  the  zonal  winds  returned  to  their 
normal  quiet-time  cycle. 

Stpler  and  Blondl  (1979)  at  Laurel  Ridge,  made  Fabry-Perot  inter¬ 
ferometer  measurements  for  the  geomagnetic  storm  of  26  Narch  1976. 
Thev  measured  equatorward  wind  of  200-S00  m/s  and  saw  a  dramatic  In¬ 
crease  in  temperature  from  750K  to  1500K.  It  is  Important  to  note  that 
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all  of  these  are  total  wind  aeasurements .  The  Beasureaents  contain 
the  normal  diurnal  wind  variation  plus  the  winds  generated  by  the 
high-latltude  source. 

Hernandez  and  Roble  (1978)  saw  evidence  of  large  scale  therao- 
spherlc  waves,  known  as  Traveling  Ionospheric  Disturbances  (TID). 
These  waves  are  atmospheric  Internal  gravity  waves  generated  by  auroral 
heating.  Roble  and  Hernandez  saw  the  meridional  winds  Increase  sharply 
in  the  equatorward  direction  with  a  large-wave  pattern  superimposed  on 
the  normal  diurnal  pattern.  The  Baxlmum  winds  were  over  300  b/s.  A 
particularly  Intense  storm  on  6  July  1974  had  equatorward  winds  of  6S0 
m/s.  Zonal  winds  also  increased  in  the  westward  direction,  reaching 
speeds  as  high  as  200  m/s.  For  the  1  April  1976  case,  by  looking  at 
the  time  difference  from  when  the  wave  was  seen  to  the  north  of  the 
station,  and  when  it  was  seen  south  of  the  station,  the  disturbance 
phase  speed  was  determined  as  400  m/s.  Four  hours  later,  a  poleward 
traveling  disturbance  with  a  phase  speed  estimated  at  630  a/s  was  seen. 
This  was  attributed  to  a  disturbance  generated  in  the  southern  hemi¬ 
sphere  auroral  zone  which  had  propagated  northward. 

There  was,  also,  an  Increase  in  temperature  from  1000K  to  1600K 
with  a  two-peaked  wave-like  structure  superimposed  on  a  mean  increase. 
Each  peak  corresponded  to  the  passage  of  a  wave  disturbance  overhead. 
The  case  on  3  May  1976  showed  considerably  more  structure  In  both  wind 
and  temperature  measurements,  indicating  a  wavelike  structure  with  a 
period  of  2  hours.  Maximum  winds  reached  400  b/s  in  the  equatorward 
direction  and  200  m/s  in  the  westward  direction. 

Roble  et  al.  (1978)  noticed  effects  of  a  gravity  wave  launched  by 


the  18  Sept  1974  sudden  commencement  on  Ionospheric  electron  density, 
electron  temperature,  and  Ion  temperature.  The  gravity  wave  passage 
caused  electron  density  to  Increase  by  as  auch  as  40Z  at  Millstone  Hill 
and  SOZ  at  Areclbo  above  qulet-tlme  levels  over  a  3-4  hour  period. 
Electron  temperature  decreased  600K  In  this  same  period.  Ion  tempera¬ 
ture  Increased  by  SOK  <SZ)  over  the  qulet-tlme  Bean. 

Mazaudler  et  al.  <1985)  using  the  Saint  Santln,  France  ISR,  showed 
that  the  thermospheric  response  to  a  geomagnetic  storm  consisted  of  an 
initial  period  of  wave  activity,  followed  by  a  longer  lasting  meri¬ 
dional  wind  disturbance  with  equatorward  winds  above  120  km,  and  low 
velocity  poleward  winds  below.  Moreover,  they  found  that  there  was  a 
time  delay  of  2-3  hours  between  the  time  geomagnetic  activity  Increased 
and  the  time  winds  above  Saint  Santln  responded  along  with  a  phase  de¬ 
lay  with  altitude.  By  subtracting  out  a  mean  qulet-tlme  wind  profile, 
they  found  that  the  equatorward  winds  due  to  the  disturbance  alone  were 
on  the  order  of  SO  m/s.  The  aagnltude,  however,  depended  on  the  level 
of  geomagnetic  activity. 

It  Is  these  aspects  of  the  thermospheric  response  to  a  geomagnetic 
storm  that  need  to  be  studied  further.  In  later  chapters,  a  generic 
storm  will  be  modelled.  This  will  be  followed  by  an  attempt  to  model 
two  of  the  storms  presented  in  Hazaudler  et  al . 


TABLE:  2.  SUMMARY  OF  STORM-TIME  OBSERVATIONS 


1.  Strong  equatorward  winds  (-300-400  m/s>. 


2.  Strong  westward  winds  (-100-200  m/s>. 


3.  Temperature  Increase  of  several  hundred  degrees  K. 


4.  Large-scale  thermospheric  waves. 


S.  A  circulation  characterized  by  a  period  of  transient  waves 
followed  bv  a  mean  circulation  with  equatorward  and  westward 
winds . 


Ill  TWO-DIMENSIONAL  THERMOSPHERIC  MODEL 


One  way  to  check  a  conceptual  model  Is  to  use  computer  simulations. 
A  computer  simulation  not  only  predicts  the  motions  of  the  thermosphere 
and  the  resultant  electrodynamics,  but  it  can  also  be  used  as  a  tool 
to  learn  about  and  better  understand  the  physics  Involved.  Addition¬ 
ally,  perturbations  introduced  into  the  simulation  can  be  used  to  de¬ 
termine  the  sensitivity  of  the  thermosphere.  This  yields  a  better  un¬ 
derstanding  of  the  Interactions  that  occur  in  the  thermosphere-iono¬ 
sphere  system. 

The  two-dimensional,  time-dependent,  thermosphere  model  of  Richmond 
and  Matsushita  (197S)  was  adopted  for  this  study.  The  model  was  ori¬ 
ginally  developed  to  simulate  the  winds  and  temperature  variations  in 
the  thermosphere  that  result  from  heating  associated  with  auroral  re¬ 
gion  electric  currents  during  a  large  hypothetical  Isolated  magnetic 
substorm.  The  model  was  later  modified  by  Blanc  and  Richmond  (1980) 
to  simulate  storm-time  ionospheric  electrical  currents. 

A  two-dimensional  model  is  a  simplification  of  reality.  Zonal 
winds  and  east-west  electric  currents  are  assumed  to  be  longitudinally 
symmetric.  The  model  is  also  symmetric  about  the  equator.  This  eli¬ 
minates  the  possibllty  of  zonal,  or  hemispheric  variations.  On  the 
otherhand,  a  two-dimensional  model  includes  much  of  the  important 
physics,  is  easier  to  setup  and  run,  and  the  results  are  easier  to  ln- 
terpet  than  a  3d  model.  (Three-dimensional  thermospheric  global  cir¬ 
culation  models  (TGCM)  have  been  developed  by  Roble  at  NCAR,  and 
Ful ler-Rowell  and  Rees  at  Univ.  College  London.  These  types  of  models 
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have  become  essential  to  a  better  understanding  of  the  thermosphere.) 


Model  Equations 


The  model  uses  the  nonllnearlzed  equations  of  thermospheric  dynam¬ 
ics  to  predict  winds,  temperatures,  pressures,  electric  fields,  and 
electric  currents.  It  is  assumed  that  the  acceleration  due  to  gravity 
<g»  is  constant  and  that  the  thermosphere  is  hydrostatic.  This  allows 
the  use  of  pressure  coordinates,  thus  eliminating  density,  which  varies 
over  6  orders  of  magnitude  in  the  thermosphere. 

The  equations  Include; 


the  hydrostatic  equation: 


+  1 

3p  p 


the  mass  conservation  equation: 


— J-r  (Vasin©)  +  =  0 
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the  zonal  component  ($)  of  the  momentum  equation: 
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the  meridional  component  (0)  of  the  momentum  equation: 
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the  energy  equation: 


3c 

3t 


1  3 


aslnd  30 


^  (v0<e  ♦  gh)sme)  +  —  (g(q  «■  V0t0  ♦  V^>  -  w(e  *  gh)) 


J0E© 


00 


and  the  equation  of  state: 


p  =  pRT  =  pgH. 

These  equations  are  fairly  standard  In  meteorological  hydrodynamics 
with  the  exception  of  the  electrodynamical  terms  such  as  JxB  in  the 
momentum  equations  and  J-E  in  the  energy  equation.  <A  list  of  symbols 
can  be  found  in  appendix  1.) 

Several  terms  incorporated  in  the  above  equations  need  to  be  ex¬ 
panded  upon.  The  energy  density  is: 

e  =  CpT  ♦  |(V|  ♦  vj>. 

The  horizontal  viscous  stress,  including  molecular  and  turbulent  co- 


TABLE:  3.  VISCOUS  AND  HEAT  CONDUCTION  PARAMETERS 


(ight 

Viscosity 

Molecular 

Turbulent 

Conduction 

Conduct io 

km 

s"1 

8_1 

s”1 

80 

1.14E-S 

2 . 89E-8 

S.67E-6 

90 

1.14E-S 

1.66E-7 

S.62E-6 

00 

1.14E-S 

7 . 44E-7 

S.43E-6 

10 

1.13E-S 

3.07E-7 

4.S6E-6 

2S 

9.64E-6 

1.14E-5 

6.93E-7 

50 

2.70E-S 

3.95E-S 

0 

!00 

1 . 00E-4 

1 . 48E-4 

0 

ISO 

3.03E-4 

S.44E-4 

0 

100 

9.S7E-4 

1 . 44E-3 

0 

ISO 

2.31E-3 

3.48E-3 

0 

100 

S.78E-3 

8.70E-3 

0 

ISO 

1.37E-2 

2.06E-2 

0 

Boundary  Conditions 


The  aodel  Is  two-dimensional  In  altitude  and  colatitude  and  forms 
a  grid  that  extends  from  the  north  pole  to  the  equator,  and  between  80 
and  450  km.  Colatitude  spacing  on  the  average  is  in  two  degree  Incre¬ 
ments  but  Is  not  constant.  It  varies  between  1.2  and  2.4  degrees,  with 
closer  spacing  In  the  auroral  zone  and  wider  spacing  at  the  pole  and 
the  equator.  The  closer  spacing  In  the  auroral  zone  allows  for  a  more 
detailed  representation  of  the  area  of  maximum  heating.  Altitude  Is 
expressed  as  a  dimensionless  height  variable  Z  =  In  <pQ/  p>  with  Z=0 
at  80  km  and  Z=14.4  at  4S0  km  for  a  total  of  49  levels  each  0.3  of  a 
scale  height  apart. 

Also,  the  horizontal  wind  velocity,  gas  energy  density,  and  pres¬ 
sure  are  held  constant  at  the  lower  boundary.  Finally,  Vq  and  go 
to  zero  at  the  pole,  and  at  the  lower  boundary.  Vq  goes  to  zero  at 
the  equator.  These  last  two  conditions  mean  that  the  boundaries  act 
like  totally  reflecting  walls. 


Initial  Conditions 


Now  that  the  basic  computational  scheme  has  been  presented,  de¬ 
tails  of  the  set-up  and  parameterization  need  to  be  discussed.  In¬ 
itially,  the  model  thermosphere  Is  at  rest.  Hence  at  subsequent  times 
the  winds,  temperatures,  and  electric  fields  and  currents  are  a  result 
of  storm  Inputs.  An  Initial  temperature  field  Is  prescribed  which 
follows  Jacchla's  1971  model  atmosphere  with  an  exospheric  temperature 
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of  1000K  <flg.  5).  Additionally,  a  background  heat  source  (Qq>  is 
specified.  This  not  an  attempt  to  model  heating  In  the  "qulet-tlme” 
thermosphere  accurately,  rather  It  provides  a  vertical  heat  distribu¬ 
tion  that  maintains  the  Initial  temperature  structure  In  the  absence 
of  auroral  heating.  Without  this  heat  source  the  model  thermosphere 
vould  slowly  cool  off. 

The  Initial  vertical  distribution  of  electron  density  Is  another 
quantity  prescribed  In  the  model.  In  this  case,  analytical  functions 
of  p  are  used  to  represent  the  electron  density  profile  In  the  E,  FI, 
and  F2  regions  based  upon  Saint  Santln  Incoherent  scatter  radar  mea¬ 
surements.  This  profile  displays  a  dual  maximum  structure  with  maxima 
at  120  km  <E  region)  and  300  km  (F  region).  The  profile  Is  constant 
with  latitude  and  represents  background  values.  (Fig. 6.) 

Electrical  conductivity  Is  calculated  using  this  profile.  (See 
chapter  V  for  a  discussion  of  conductivity. )  The  Pedersen  conductivity 
(  )  and  the  Hall  conductivity  (  >  are: 


0^(2, 0,t>  =  o^(z,0,t>  r(z,0> 

where  N(z,0,t>  Is  the  electron  density,  e  is  the  electron  charge,  and 
r(z,0>  Is  the  collision  frequency  divided  by  the  electron  gyrofre- 
quency.  B(0)  is  the  dipolar  magnetic  field  which  Is  aligned  with  the 


earth's  rotational  axis. 
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and  are  Important  In  the  calculation  of  currents  and  total  power. 
Values  of  o^  and  o^  are  dependent  on  both  the  background  electron  den¬ 
sity  and  enhancements  due  to  particle  precipitation  during  a  storm. 

The  driving  force  In  the  model  is  the  Input  of  energy  Into  the  au¬ 
roral  zone.  There  are  three  major  mechanisms  that  can  provide  this 
energy:  Ion  drag  forces,  particle  precipitation,  and  Joule  heating. 
The  details  of  these  mechanisms  were  presented  earlier.  The  goal  here 
Is  to  discuss  how  they  are  parameterized  In  the  model. 

Ion  drag  Is  only  weakly  coupled  to  the  storm  time  winds  In  the 
model .  The  appropriate  JxB  terms  are  incorporated  In  the  momentum 
equations,  but  the  electrical  conductivity  used  In  the  calculation  of 
J  doesn't  Increase  during  the  storm.  This  Is  a  result  of  not  allowing 
electron  density  to  Increase,  and  was  done  because  of  the  zonal  sym¬ 
metry  of  the  model.  Typically,  during  a  large  storm,  such  as  the  one 
modelled  here,  the  effects  of  Ion  drag  would  switch  direction  with  time 
as  the  earth  rotates  underneath  the  magnetospherlc  convection  pattern. 
In  a  zonally  symmetric  model  there  Is  no  direct  way  to  represent  local 
time  dependences,  so  the  model  provides  a  background  value  for  Ion  drag 
and  Ignores  the  storm-time  Influence. 

The  model  parameterizes  particle  precipitation  by  an  analytical 
function  but  only  as  It  pertains  to  electron  density  for  Joule  heat¬ 
ing.  Hence,  the  effect  of  an  enhanced  electron  density  on  the  overall 


conductivity,  or  heating  due  to  the  particles  themselves  Is  Ignored. 


Blanc  and  Richmond  <1980)  estimated  thst  particle  precipitation  was 
small  compared  to  Joule  heating.  Also,  the  energetic  particles  pene¬ 
trate  through  the  thermosphere  to  the  E  region  where  the  density  Is 
high  and  air  motions  are  difficult  to  Initiate. 

Joule  heating  Is  the  most  Important  heat  source  In  the  model.  How 
well  the  model  parameterizes  Joule  heating  Is  therefore  crucial.  The 
first  question  to  answer  Is  where  does  the  heating  occur?  In  the 
thermosphere  the  maximum  heating  takes  place  In  the  auroral  zone.  The 
model  parameterizes  this  spatial  characteristic  by  using  a  guasslan 
that  peaks  at  21°  colatitude  and  tails  off  to  zero  at  the  edges  of  the 
auroral  zone  (Figure  7).  (The  following  two  paragraphs  are  from  work 
done  with  C.  Hazaudler  and  A.D.  Richmond,  Mazaudler  et  al.,1987.) 

The  second  question  to  answer  Is  how  to  parameterize  the  magnitude 
of  the  Joule  heating?  One  way  is  to  use  Incoherent  scatter  radar  ob¬ 
servations.  Incoherent  scatter  radar  measures  ion  density  and  Ion 
velocltv.  From  these  quantities,  electric  fields  and  conductivities 

can  be  derived.  A  direct  estimation  of  Joule  heating  can  be  made  us- 

2 

lng  the  expression  U.  =  a  E  ,  where  a  Is  the  height  Integrated 

J  P  P 

Pedersen  conductivity  and  E  Is  the  perpendicular  electric  field.  Brekke 

Mq7B>,  Banks  (1977),  Banks  et  al .  (1981),  Vickrey  et  al.  (1982), 

Duboin  and  Kamlde  (1984)  provide  estimations  of  Joule  heating  varying 
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from  2  to  50  ergs  cm  s  .  The  problem  with  using  Incoherent  scatter 
radar  observations  Is  that  they  only  provide  Joule  heating  rates  over 
a  small  area,  end  only  during  the  daytime. 

A  second  method  of  parameterizing  Joule  heating  Is  by  using  mag- 
netir  data.  Global  Ionospheric  electric  fields  and  electric  current 


distributions  can  be  determined  through  Inversion  of  magnetic  data 


(Kamlde  et  al.,  1981;  Karalde  et  al.,  1982).  The  electric  field  and  a 
conductivity  model  (Kamlde  and  Richmond,  1982)  determine  Joule  heating. 
This  estimation  of  the  Joule  heating  Is  not  as  accurate  as  the  estima¬ 
tion  from  Incoherent  scatter  radar  measurements,  but  It  does  provide 
a  global  estimate.  Comparisons  between  incoherent  scatter  estimates 

of  Joule  heating  (U^ )  and  variation  of  the  H  component  of  the  magnetic 
2 

field  (AH)  at  the  same  location  have  been  made  by  Brekke  and  Rlno 

(1978)  and  Duboln  and  Kamlde  (1984)  in  order  to  define  a  coefficient 

2 

h'  such  that  =  h'(AH)  These  papers  show  that  the  magnitude  of 

the  h'  parameter  changes  with  the  chosen  magnetic  event.  This  under¬ 
scores  the  point  that  Joule  heating  Is  not  simply  related  to  the  mag¬ 
netic  field  variation.  Ahn  et  al.  (1983),  and  Baumjohann  and  Kamlde 
(1984)  made  comparisons  between  the  AE  index  and  Joule  heating  derived 
from  magnetic  data,  In  order  to  define  the  Joule  heating  rate  Integrat¬ 
ed  over  one  hemisphere  as  a  function  of  the  AE  index.  Table  4.  shows 
the  different  expressions  found  by  these  authors.  The  table  shows  two 


formulations  of  this  expression  from  Ahn  et  al.  based  on  AE  derived 


from  71  stations  In  the  first  case,  and  from  12  stations  In  the  second 
case.  The  expression  from  Baumjohann  and  Kamlde  is  based  on  71  sta¬ 
tions.  The  main  difference  between  these  papers  is  the  conductivity 
model  that  was  used.  In  this  study  a  hybrid  expression  Is  used,  based 
on  Baumjohann  and  Kamlde  with  the  factor  multiplying  the  AE  Index  In¬ 
creased  by  0.2  to  account  for  the  use  of  AE(12)  rather  than  AE(71). 
(In  accordance  with  Ahn  et  al.)  This  yields  a  hemispheric  heating  rate 
of  U  .  < 12 >  =  4.0  x  108  AE< 12 >  Watts. 


Ahn  etaL  1983 


Uj  =  2.3  x  10*  AE(12)  in  Watts 


Ahn  et  aL  1983  Uj  =  1.9  x  10®  AE(71)  in  Watts 

Baumjohann  and  Kamide  Uj  =  3.3  x  10s  AE(71)  in  Watts 
1984 


Used  in  Mode!  (Hybrid)  Uj  =  4.0  x  10*  AE(12)  in  Watts 

AE(12)  -»  index  derived  from  12  stations 
AE(71)  ^  index  derived  from  71  stations 


Table  4.  Expressions  relating  hemispheric  Joule  heating  to  the  AE 
Index.  AE<12>  and  AE(71)  are  the  Auroral  Electrojet  Index  in  nano¬ 
teslas  based  on  12  and  71  stations  respectively. 


Computational  Technique 


Two  types  of  computations  are  made  In  the  model.  They  are  the  cal¬ 
culation  of  flux  divergences  and  the  time  Integration  of  the  dependent 
variables.  Flux  divergences  are  calculated  from  finite  differences, 
eg: 
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where  x  Is  a  flux  and  x^+^  and  are  values  of  x  to  either  side  of 

where  the  flux  Is  being  calculated.  Ap^  Is  the  thickness  of  the  1th 
pressure  layer. 

Three  quantities  are  stepped  foward  in  time.  They  Include  the 

meridional  wind  (V  ),  the  zonal  wind  (V  )  and  the  energy  flux  density 
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(e),  which  Is  used  to  predict  temperature.  The  time  integration  is 
done  as  follows  (Richmond  and  Matsushita,  1975): 

The  momentum  equations  for  V  or  V  or  the  energy  equation  can  be 
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and  solved  for  V(t+At)  by  Iteration  with  H*  being  the  1th  Iterated 
value  of  V<t+At>.  The  zero  Iteration  Is: 


V  <t+At  >  =  2  V(t>-  V(t-At). 


Then: 


V1  ( t+At )  =  H»(t)  ♦  At*F[  ~  V(t)  ♦  ~  «p1-1<t+At>,  t  ♦  |  At]  J  . 


Richmond  and  Matsushita  <197S>  found  that  only  2  Iterations  were 
needed  to  converge  on  an  answer  and  that  a  time  step  of  30sec  was 
sufficient  to  maintain  computational  stability.  For  each  time  step 
this  method  is  used  for  each  of  the  quantities  at  each  element  of  the 
grid. 


IV  MODEL  RESULTS 


In  order  to  simulate  the  thermospheric  response  to  a  geomagnetic 
storm,  an  energy  Input  profile  was  devised  to  provided  an  Input  for 
the  two-dimensional  model.  The  profile,  which  is  shown  in  figure  8., 
represents  a  generic  storm.  The  abrupt  swltch-on  of  energy  is  similar 
to  what  would  be  expected  from  a  sudden  storm  commencement.  This  re¬ 
sults  In  an  initial  period  of  transient  waves,  followed  by  the  develop¬ 
ment  of  a  mean  circulation.  "he  storm  is  maintained  at  a  constant 
power  level  of  4X1011  Watts,  for  a  period  of  12  hours.  This  Is  rough¬ 
ly  equivalent  to  a  geomagnetic  storm  with  an  AE  of  1000  for  the  same 
period. 

The  power  level  of  the  generic  storm  represents  an  average  storm. 
It  Is  on  the  same  order  of  magnitude  as  the  estimates  of  Joule  heating 
for  the  three  storms  presented  by  Hernandez  and  Roble  (1978)  and  the 
sudden  storm  commencement  discussed  by  Roble  et  al .  (1978).  It  falls 
between  values  of  Joule  heating  for  the  "non-storm"  of  Dickinson  et 
al.  <1975)  (0. 6-1 . OXIO^W)  and  a  very  strong  storm  (8.24X10**W)  used 
by  Blanc  and  Richmond  (1980). 

What  Is  unusual  about  this  profile.  Is  the  period  for  which  this 
storm  is  maintained.  Typically,  observed  storms  (and  most  of  the 
modelled  storms)  have  enhancements  in  Joule  heating  that  last  for  3-4 
hours.  Hence,  12  hours  Is  a  long  time  to  maintain  this  level  of  heat¬ 
ing.  Also,  a  typical  storm  has  a  heating  rate  that  fluctuates  with  a 
period  of  2-4  hours.  By  maintaining  the  Joule  heating  at  this  rate 
for  12  hours,  the  storm-time  mean  circulation  Is  able  to  approach  a 
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steady-state  (although,  It  still  hasn't  been  reached  at  this  point). 
The  addition  of  periodicity  In  the  heating  profile  produces  a  similar 
periodicity  In  the  response  (eg.  a  forced  oscillator).  Although,  as 
Francis  (1975)  has  stated  the  oscillatory  nature  of  large  traveling 
disturbances  does  not  necessarily  Imply  an  oscillatory  source.  The 
simple  heating  profile  used  here,  however,  makes  It  easier  to  separate 
out  the  effects  of  the  transient  waves  and  the  mean  circulation. 

Model  Generated  Waves 

The  Impulsive  heating,  which  results  from  this  heating  profile, 
initiates  strong  upward  motion  that  generates  waves.  These  waves  show 
up  as  a  perturbation  In  the  meridional  wind  field  that  starts  at  the 
heat  source  In  the  auroral  zone  and  propagates  both  poleward  and  equa- 
torward.  Figure  9.  shows  the  propagation  of  the  disturbance  for  vari¬ 
ous  times  at  an  altitude  of  250  km.  This  altitude  was  chosen  because 
it  was  in  the  middle  of  the  model  grid  and  was  the  same  altitude  as  the 
majority  of  the  observations. 

The  wind  disturbance  arrives  at  the  equator  in  3-4  hours  with  an 
amplitude  of  55  m/s  in  the  equatorward  direction.  This  is  similar  to 
the  amplitudes  observed  by  Mazaudler  et  al.  (1985).  However,  even 
after  subtracting  out  a  quiet-time  mean  wind,  the  winds  seen  by 
Hernandez  and  Roble  (1978)  are  still  larger  by  a  factor  of  three. 
Since  the  model  gives  a  linear  response  to  heat  inputs,  it  would  take 
a  storm  three  times  as  large  as  used  here  to  achieve  a  similar 
amplitude.  The  disturbance  also  exhibits  a  phase  speed  of  524  m/s 


(§250  km),  which  is  in  line  with  that  observed  by  Hernandez  and  Roble 
<1978) . 

The  disturbance  is  actually  a  wave  packet  made  up  of  a  whole  spec¬ 
trum  of  waves  (Hines,  1974)  of  various  periods  and  wavelengths.  This 
makes  It  difficult  to  determine  a  specific  period  or  wavelength.  The 
model -produced  waves  exhibit  a  period  on  the  order  of  100  minutes  and 
a  wavelength  between  2000  and  4000  km.  The  packet  appears  as  a  single 
Impulse  at  the  head  of  the  disturbance  because  the  waves  that  form  the 
oscillatory  tail  of  the  packet  have  been  dissipated  via  thermal  conduc¬ 
tion  and  viscosity.  (Chlmonas  and  Hines,  1970;  and  Richmond  and 
Matsushita,  197S.)  The  waves  appear  to  be  Increasing  in  period  at 
lower  altitudes  and  with  Increasing  distance  from  the  source  region. 
This  is  a  consequence  of  the  fact  that  the  angle  at  which  the  wave  is 
preferentially  launched  is  dependent  on  frequency.  The  longer  the 
period  of  the  wave,  the  smaller  this  angle  is.  Hence,  the  longer 
period  waves  would  be  predominantly  at  lower  altitudes  and  distances 
farther  away  from  the  source  (Richmond  and  Matsushita,  1975;  Francis, 
1975 ) . 

The  equatorward  propagation  of  the  waves  can  also  be  seen  in  terms 
of  a  temperature  perturbation.  Figure  10.  shows  similar  wave  features 
to  those  seen  in  the  meridional  winds,  but  in  this  case  the  amplitude 
represents  a  4Z  increase  in  temperature.  This  is  an  increase  of  38K 
over  an  initial  temperature  of  942K.  Roble  et  al.  (1978)  showed  a 
similar  perturbation  in  neutral  temperature.  Also,  after  the  wave  has 
passed,  the  temperature  has  Increased  by  35K  at  the  equator,  signifying 
that  the  waves  have  deposited  energy  at  low  latitudes. 


Figure  9.  Storm  generated  meridional  winds  at  250  Km.  Storm  time 
<0-5  hrs)  vs  colatitude  is  shown.  Each  tic  mark  represents  25  m/s 
positive  southward. 


gure  10.  Same  an  Figure  9.  except  for  temperature  change.  Fach 
c  mark  represents  1Z  temperature  change  from  the  initial  background 
mperature  <942K>. 


Figures  11.  and  12.  show  the  disturbance  as  It  would  be  seen  at  a 
single  colatitude  for  the  Meridional  velocity  and  temperature  change, 
respectively.  This  16  how  the  observing  stations  would  see  the  dis¬ 
turbance  A  colatitude  of  48  degrees  was  chosen  because  it  corresponds 
to  the  colatitude  at  which  the  majority  of  the  observations  are  made. 
Above  2S0  km,  the  waves  are  pretty  much  in  phase  with  altitude.  This 
is  due  to  the  dominant  Influence  of  viscosity  on  the  bulk  motion  of  the 
air  Below  2S0  km,  the  waves  exhibit  a  distinct  phase  lag  with  alti¬ 
tude;  this  lag  Is  about  1.66  hours  from  2S0  km  to  120  km.  This  phase 
lag  leads  one  to  believe  that  the  waves  are  gravity  waves.  The  phase 
tilt  Is  what  Is  expected  of  gravity  waves  with  a  low  altitude  source 
(120  km)  and  with  energy  (group  velocity)  propagating  upward  (Hines, 
1960).  Figure  12.,  also  shows  that  after  five  hours,  there  Is  a  gen¬ 
eral  increase  in  temperature,  resulting  from  the  development  of  the 
mean  circulation. 

Mean  Circulation 

Heating  In  the  auroral  zone  leads  to  the  development  of  a  pressure 
bulge  at  high  latitudes  and  hence  a  pressure  gradient  that  drives  hori¬ 
zontal  winds.  Initially,  there  Is  a  period  of  transient  gravity  waves 
as  described  In  the  previous  section.  After  these  waves  have  passed 
through  the  system,  a  mean  circulation  developes. 

Figures  13a, b,  and  c  show  the  model  outputs  for  the  meridional, 
zonal,  and  vertical  components  of  the  thermospheric  wind,  respectively, 
two  hours  Into  the  simulation.  Figure  13d,  shows  temperature  changes, 
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also  after  two  hours.  <  It  is  Important  to  note  that  the  transient 
waves  are  evident  at  two  hours;  however,  the  beginnings  of  the  mean 
circulation  can  also  be  seen  at  this  time.) 

Figure  13a  shows  that  the  leading  edge  of  the  wind  disturbance  has 
reached  70°  colatitude  at  the  upper  boundary  of  the  model.  At  lower 
altitudes,  the  disturbance  lags  behind  this.  Just  equatorward  of  the 
auroral  zone,  the  developing  circulation  Is  confined  between  the  alti¬ 
tudes  of  120  and  200  km,  and  between  colatitudes  of  22.5°  and  4S°  with 
maximum  equatorward  winds  of  110  m/s.  Poleward  of  the  auroral  zone, 
there  is  a  smaller  poleward  wind  with  maximum  winds  of  30  m/s.  Above 
200  km,  the  transient  gravity  waves  (as  described  earlier)  are  domi¬ 
nant  . 

Figure  13b,  shows  the  zonal  winds  at  the  two  hour  point.  These 
winds  are  the  result  of  the  corlolis  force  acting  on  the  meridional 
winds.  Ion  drag  forces,  as  shown  by  Richmond  and  Matsushita  (1975), 
would  also  affect  the  zonal  component  of  the  winds.  However,  ion  drag 
forces  are  suppressed  in  this  version  of  the  model  for  reasons  discuss¬ 
ed  earlier.  Westward  zonal  winds  extend  from  22. S°  to  45°  colatitude 
above  120  km  with  a  maximum  of  50  m/s.  Poleward  of  the  auroral  zone, 
there  are  eastward  zonal  winds,  which  reach  a  maximum  of  40  m/s. 

The  vertical  component  of  the  wind  (Figure  13c)  at  this  time  shows 
a  core  of  upward  moving  air  centered  about  21°  colatitude  with  maximum 
vertical  velocities  of  8  m/s.  Equatorward  of  this  core  to  76°  colati¬ 
tude,  the  vertical  motion  is  downward  with  velocities  less  than  1  m/s, 
except  above  200  km  where  gravity  wave  activity  is  present. 

Figure  13d  shows  a  temperature  increase  over  the  model's  initial 


temperature  field.  There  has  been  a  significant  temperature  Increase 
in  the  region  between  the  pole  and  40°  colatitude,  and  between  110  and 
170  km.  The  maximum  Increase  Is  110K  at  135  km  between  the  pole  and 
the  auroral  zone.  At  40°  colatitude,  the  temperature  has  increased  40K 
at  135  km.  The  temperature  increase  is  the  consequence  of  Joule  heat¬ 
ing  at  E  region  heights  in  the  auroral  zone,  which  the  winds  then 
transport  away  from  the  source.  Above  200  km,  and  centered  about  the 
auroral  zone  Is  another  region  of  Increasing  temperature.  This  in¬ 
crease  Is  the  result  of  Joule  heating  at  F  region  heights  and  has 
reached  70K  after  two  hours.  Above  200  km,  temperature  perturbations 
are  associated  with  the  transient  gravity  waves  are  seen.  Figure  13d, 
also  shows  the  equatorward  extent  of  the  disturbance  (colatitude  70° 
at  the  model's  upper  boundary). 

Figure  14a  shows  that  after  5  hours,  the  effects  of  the  mean  cir¬ 
culation,  centered  about  150  km,  have  reached  the  equator.  Just  equa¬ 
torward  of  the  auroral  zone  (25°  colatitude),  the  maximum  equatorward 
winds  are  130  m/s.  At  mid-latitudes,  the  winds  are  equatorward  at  50- 
70  m/s,  which  is  in  reasonable  agreement  with  Mazaudier  et  al.  (1986) 
and  other  storm-time  observations  (once  a  quiet-time  mean  is  sub¬ 
tracted).  At  low  latitudes,  the  equatorward  winds  are  10-30  m/s. 
Notice  also,  that  the  winds  seem  to  be  confined  to  heights  above  120 
Km.  This  is  due  to  the  higher  air  density  below  120  km.  In  fact,  a 
weak  poleward  wind  developeB  (less  than  10  m/s)  at  the  lower  altitudes 
to  satisfy  the  need  to  conserve  mass  in  the  system  (a  return  flow). 

Above  240  km,  the  effects  of  the  waves  have  weakened,  but  are  still 


evident.  At  5  hours,  there  appears  to  be  a  poleward  propagating  wave 
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at  low  latittides.  This  wave  1b  an  artifact  of  the  model  boundary  con¬ 
ditions.  Since  V  must  go  to  zero  at  the  equator,  the  equatorward 

♦ 

boundary  acts  like  a  reflecting  wall.  Hence,  waves  that  were  origin¬ 
ally  propagating  toward  the  equator  are  reflected  back  towards  the 
pole.  Some  physical  significance  can  be  attached  to  this  reflected 
wave.  It  represents  a  wave  that  was  generated  In  the  southern  hemi¬ 
sphere's  auroral  zone  and  has  propagated  northward  Into  the  northern 
hemisphere . 

After  S  hours,  the  westward  winds  have  expanded  equatorward  to  76° 
colatitude  and  the  maximum  speed,  located  at  ISO  km  and  34°  colatitude, 
has  Increased  to  120  m/s  (figure  14b>.  The  vertical  velocities  (figure 
14c >  show  that  the  wave  activity  has  all  but  disappeared,  leaving  a 
"Jet-llke"  core  of  upward  vertical  motion  centered  at  21°  colatitude. 
This  core  exhibits  a  double-maxima  structure,  one  In  the  E  region  (with 
a  max  wind  of  4  m/s)  and  one  in  the  F  region  (with  a  max  wind  wind  of  8 
m/s).  Elsewhere,  the  vertical  winds  are  predominately  downward  with 
velocities  less  than  1  m/s.  Also,  at  this  time,  the  temperature  at 
low  latitudes  has  Increased  by  20-60K.  In  the  auroral  region  at  130 
Km,  the  temperature  has  Increased  by  140K  (figure  14d). 

Three  hours  later,  after  the  storm  has  been  going  for  8  hours,  the 
meridional  winds  have  maintained  their  maximum  speed  but  their  spatial 


extent  seems  to  have  decreased  (figure  15a).  This  may  be  due  to  the 
fact  that  the  Initial  temperature  gradient  (which  Is  driving  the  winds) 
has  been  changed  by  the  winds  themselves  and  the  circulation  is  now 
readjusting.  The  zonal  winds  (figure  ISb)  have  continued  to  increase 


In  magnitude  with  a  maximum  of  140  m/s.  Their  horizontal  extent  has 
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also  Increased  equatorward  to  77°  colatitude.  The  vertical  velocities 
have  Increased  slightly  In  magnitude  (figure  15c),  and  the  temperature 
at  high  latitudes  has  Increased  by  as  much  as  160K  while  at  low  lati¬ 
tudes  It  has  Increased  by  60K  (figure  lSd). 

After  11  hours,  the  mean  circulation  approaches  a  steady-state 
(although,  it  hasn't  been  reached  yet).  The  meridional  winds  (figure 
16a)  has  changed  little  over  the  last  three  hours,  except  they  have 
expanded  in  spatial  extent  slightly.  The  zonal  and  vertical  winds 
have  changed  little  over  the  last  three  hours  (figure  16b, 16c).  A 
small  area  of  eastward  winds  has  developed  equatorward  of  the  auroral 
zone  and  below  120  km.  This  may  be  associated  with  the  coriolis  force 
acting  on  the  low  altitude  poleward  winds  (the  return  flow). 

Finally,  figure  16d  shows  that  the  thermospheric  temperature  in  the 
auroral  zone  has  increased  by  as  much  as  200K.  At  low  latitudes  tem¬ 
peratures  have  increased  by  as  much  as  100K  due  to  dynamically  tran¬ 


sported  heat  from  the  auroral  zone  to  the  equator. 


Figure  13a.  Meridional  winds  (positive  equatorvard)  2  hours  Into 
storm  simulation.  Contours  are  In  10  m/s  Increments. 


lgnre  13b.  Zonal  winds  (positive  eastward)  Z  hours  into  storm 
limitation.  Contours  are  In  10  m/s  Increments. 


Figure  13c.  Vertical  velocities  (positive  upward)  2  hours  Into  storm 
simulation.  Contours  are  in  1  m/s  Increments. 
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Figure  Ibb.  Same  as  Figure  13b.,  except  at  11  hours 


13d.,  except  at  1 1  hours 


V  ELECTRODYNAMICS 


Thermospheric  winds  move  charged  particles  in  the  ionosphere 
through  the  earth's  magnetic  field.  This  results  in  the  generation  of 
electromotive  forces  which  drive  electrical  currents,  and  set  up  elec¬ 
tric  fields  (Richmond,  1979a). 

In  general,  moving  charged  particles  are  acted  upon  by  the  mag¬ 
netic  Lorentz  force  <-q  $  x  i  )  and  gyrate  around  magnetic  field  lines. 
This  restricts  the  movement  of  particles  in  the  direction  perpendicular 
to  field  lines.  Motion  along  the  field  lines,  however.  Is  not  acted 
upon  by  magnetic  forces  and  charged  particles  therefore  can  easily  move 
along  them. 

A  force  is  required  to  move  charged  particles  across  field  lines. 

A  perpendicular  electric  field  imposes  such  a  force.  At  altitudes 
above  180  km,  both  electrons  and  ions  move  together  with  the  f  X  ? 

drift.  Since  there  is  no  charge  separation,  there  is  no  current.  Be¬ 
low  80  km,  frequent  collisions  between  charged  particles  and  neutrals 
cause  the  charged  particles  to  move  together  with  the  neutral  wind, 
and  again  there  is  no  current.  Above  80  km,  electrons,  because  of 
their  small  size  and  mass,  undergo  few  collisions  and  move  with  the 
?  X  ^  drift.  Ions,  on  the  otherhand,  have  frequent  collisions,  which 
cause  their  motion  to  deviate  from  the  2  X  £  direction.  Collisions 
affect  ion  motion  in  two  ways.  First,  they  stop  an  ion's  foward 
motion,  which  allows  the  electric  force  to  accelerate  the  ion  in  the 
direction  of  the  electric  field.  The  ion  now  has  components  of  motion 
in  the  2  X  2  direction  and  in  the  direction  of  2.  This  causes  charge 


separation  (as  the  electrons  move  only  In  the  2  X  S  direction)  and  a 
Pedersen  current  in  the  direction  of  2  flows.  Second,  collisions  re¬ 
tard  the  Ion's  motion  In  the  £  X  $  direction  relative  to  the  electrons, 
causing  charge  separation  and  a  Hall  current  in  the  X  5  direction 
(the  electrons  "carry"  the  current). 

The  currents  can  be  summerlzed  by  ionospheric  Ohm's  law: 

3  =  aQ  2  ♦  ot  (  ^  X  3  )  +  a2  b  X  (  ^  +  lx  I  ) 

The  first  term  Is  the  parallel  current,  where  a Is  the  parallel  con¬ 
ductivity.  The  second  term  Is  the  Pedersen  current,  where  is  the 

Pedersen  conductivity.  The  third  term  Is  the  Hall  current,  where  a 

* 

is  the  Hall  conductivity  and  b  Is  the  unit  vector  In  the  direction  of 

a. 

Conductivity  relates  to  the  mobility  of  charged  particles,  and  Is 
dependent  on  the  ratio  of  the  collision  frequency  (\>)  to  the  gyrofre- 
quency  (u»).  The  parallel  conductivity  Is  : 

°0  =  ^e  +  ^1  I  e  I  ^ 

®eue  ®lul  . 

The  parallel  conductivity  increases  with  altitude  and  quickly  becomes 
large.  This  allows  currents  to  flow  along  field  lines,  which  shorts 
out  2 1|  Hence,  this  term  can  be  ignored  In  Ionospheric  Ohm’s  law. 
Additionally,  the  field  lines  act  like  equlpotenttals  and  2  Imposed  by 
the  magnetosphere  can  be  mapped  down  to  the  ionosphere. 

The  Pedersen  and  Hall  conductivities  are: 
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and 


respectively.  The  Pedersen  conductivity  maximizes  at  140  km,  where 
<i>j  =  vj  j  and  is  In  the  direction  of  2.  The  Hall  conductivity  max¬ 
imizes  at  10S  km,  where  =  Ueue>  an<*  18  *n  t^ie  ^  X  £  direction. 
Pedersen  and  Hall  currents  then  arise  from  either  an  electrostatic 
field  or  an  Induced  electric  field,  as  described  by  ionospheric  Ohm's 
law. 

The  qulet-tlme  diurnal  and  semi-diurnal  tidal  motions  of  the 
thermosphere  drive  electrical  currents.  These  currents  make  up  the  Sq 
current  system,  and  can  be  determined  from  ground-level  magnetic 
varlations< f ig  17).  On  the  dayslde,  the  result  Is  two  large  current 
vortices.  In  the  northern  hemisphere,  the  current  flows  counterclock¬ 
wise,  and  In  the  southern  hemisphere  the  current  flows  clockwise.  The 
maximum  current  Is  on  the  order  of  6  X  104  Amps.  On  the  nightside, 
where  the  conductivity  is  lower,  the  currents  are  half  the  strength  of 
the  dayslde  currents,  and  flow  in  the  opposite  direction. 

Disturbance  Dynamo  Model 

Storm-time  thermospheric  winds  can  also  drive  electrical  currents. 


Auroral  heating  due  to  geomagnetic  storms  sets  up  a  large-scale 
"Hadley-cell"  like  circulation  with  winds  flowing  southward  at  mld- 
latltudes  above  120  km  (as  described  In  Chapter  IV).  Conservation  of 
angular  momentum  results  In  the  development  of  a  westward  component  In 
the  wind.  This  westward  wind  seta  up  an  Induced  electric  field  ( 
x  i,  t  downward  In  the  northern  hemisphere)  that  yields  an  equatorvard 
Pedersen  current.  This  leads  to  the  buildup  of  positive  charges  at  the 
equator  and  negative  charges  at  high  latitudes,  which  results  In  a 
poleward  polarization  electric  field  Just  large  enough  to  balance  the 
equatorvard  current.  The  poleward  electric  field  interacts  with  the 
earth’s  magnetic  field  (&  down)  leading  to  an  eastward  Hall  current. 
When  day/night  conductivities  are  considered,  the  eastward  Hall  current 
Is  stopped  at  the  terminator.  Closure  of  the  current  leads  to  a 
current  vortex  similar  to  the  Sq  current  system  (on  the  same  order  of 
magnitude),  but  with  the  opposite  sense.  This  is  the  antl-Sq  current 
system  (Blanc  and  Richmond,  1980).  (Figure  18.) 

The  model  generates  both  the  northward  electric  field  and  the 
eastward  current.  Figure  19.  shows  the  development  of  the  electric 
field.  Once  the  storm  starts,  the  electric  field  starts  growing  and 
reaches  a  maximum  of  4.S2  «V/m,  shortly  after  12  hours.  After  storm 
shutoff,  the  electric  field  slowly  starts  to  decrease  in  strength. 
The  slow  decrease  results  from  the  slow  dissipation  of  the  zonal  winds, 
which  continue  to  drive  currents  and  electric  fields  even  after  the 
storm  has  stopped.  (This  was  referred  to  by  Lyons  et  al .  (198S)  as 
the  neutral  wind  "flywheel".)  Figure  20.  shows  the  eastward  Hall  cur¬ 
rent.  The  current  also  builds  up  in  strength  as  the  storm  progresses, 


reaching  a  maximum  of  S6.S  mA/m  after  12  hours  (approx.  1.25x10  Amps, 
which  Is  1/5  the  magnitude  of  the  Sq  current  at  44°  colatitude).  Like 
the  electric  field,  the  current  slowly  decreases  after  the  storm  ends. 
In  the  next  chapter,  the  model  produced  electric  field  will  be  compar¬ 
ed  to  measurements. 
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VI  COMPARISON  BETWEEN  SAINT -SANT IN  OBSERVATIONS  AND  MODEL 
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PREDICTIONS 

In  this  chapter,  neutral  wind  and  electric  field  data  fro*  two 
storms  analysed  by  Mazaudler  et  al.  (198S)  will  be  compared  to  the 
model.  (This  section  Is  from  work  done  with  C.  Mazaudler  and  A.D. 
Richmond,  Mazaudler  et  al.,1987.)  The  two  storms  are  June  2,  1978, 
and  April  19,  1977  (figures  21  and  23).  The  first  storm  shows  good 
agreement  between  the  observations  and  the  simulations,  while  the  se¬ 
cond  storm  shows  only  fair  agreement.  The  top  panels  In  figures  21. 
and  23.  show  the  AE  Index  for  these  two  storms.  AE  represents  the 
time  history  of  Joule  heating  used  in  the  simulations.  The  middle  and 
bottom  panels  show  the  storm  Induced  southward  neutral  winds  derived 
from  observations  and  the  model  produced  winds.  Note  that  the  alti¬ 
tudes  used  In  the  first  case  differ  from  the  second  case.  Figure  24. 
compares  the  observed  electric  field  with  model  simulations  for  the 
April  19  case.  (The  electric  field  was  unavailable  for  the  June  1978 
storm  due  to  technical  problems  with  the  radar.) 

The  Saint  Santin  Incoherent  scatter  radar  was  discussed  In  an 
earlier  chapter.  Not  mentioned  was  the  fact  that  the  time  resolution 
of  the  wind  measurements  is  30  minutes  with  an  accuracy  of  1  15  m/s. 
In  addition,  it  Is  important  to  remember  that  the  observed  storm- 
induced  winds  are  obtained  by  subtracting  the  quiet-time  winds 
(Mazaudler  et  al.,  1985)  from  the  observations.  The  observed  winds 
exhibit  much  day  to  day  varlablity  even  on  magnetically  quiet  days. 
This  varlablity  Is  not  well  understood  and  complicates  the  interpeta- 
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tlon  of  storm-time  variations. 

Case  #1  June  2,  1978 

The  middle  panel  of  figure  21.  shows  a  comparison  between  the  neu¬ 
tral  winds  derived  from  observations  and  the  neutral  winds  predicted 
bv  the  model  at  150  km.  The  observations  show  several  peaks  on  the 
order  of  70  to  80  m/s.  The  first  peak  from  9-10UT  corresponds  to 
the  first  peak  In  AE.  The  second  peak  In  the  neutral  wind  between  12 
and  14UT  corresponds  to  the  second  peak  In  the  AE  Index.  Then  beginn¬ 
ing  at  17-18UT,  there  Is  a  possibility  of  a  third  peak  In  the  neutral 
winds  before  the  observations  terminate  at  dusk.  The  AE  Index  shows  a 
large  Increase  from  14-20UT,  but  one  can  only  speculate  at  what  Is 
happening  to  the  neutral  wind.  This  case  shows  a  3-4  hour  time  lag 
between  Increases  In  the  AE  Index  and  Increases  In  the  neutral  wind. 

For  this  case,  the  computer  model  shows  good  qualitative  agreement 
with  the  observations.  The  simulated  winds  Bhow  a  southward  wind 
component  of  about  the  same  magnitude  as  the  observations,  although 
the  time  variations  are  different.  The  time  lag  between  peaks  In  the 
model  neutral  winds  and  the  AE  index  Is  about  one  hour  shorter  than 
the  time  lag  between  the  observations  and  the  AE  Index.  The  reason 
for  this  difference  Is  not  understood 

Figure  22.  shows  a  vertical  profile  of  the  meridional  winds.  On 
the  lefthand  side  is  the  qulet-tlme  average  for  June  1-14.  On  the 
rlghthand  side  Is  a  comparison  between  the  model  wind  profile  and  the 
observations  averaged  for  June  2.  Above  120  km,  the  model  and  the 
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Figure  22.  Vertical  profiles  of  daily  mean  Saint  Santin  southward 
neutral  winds.  Lefthand  side  shows  "quiet-time  average"  for  June 
1-14,  1978.  Rlghthand  side  shows  a  comparison  between  the  obser¬ 
vations  and  the  model  for  June  2,  1978. 
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observations  show  excellent  agreement.  Below  120  km,  the  model  shows 


a  slight  northward  flow,  while  the  observations  show  a  substantial 
northward  flow  down  to  10S  km.  In  the  model,  this  northward  flow  re¬ 
presents  the  return  mass  flow,  which  balances  the  southward  flow  at 
higher  altitudes.  The  observed  wind  Is  probably  not  perfectly  sym¬ 
metric  In  longitude,  hence  the  average  meridional  mass  flow  Is  not  re¬ 
quired  to  balance  out  between  different  heights,  as  It  must  In  the 
model.  This  allows  for  the  observed  winds  at  lower  altitudes  to  differ 
from  the  model  winds.  Also,  the  observed  winds  below  120  km  could  be 
the  result  of  the  normal  day  to  day  varlabllty  in  the  winds  and  not 
storm-related  effects. 

Case  #2  April  19,  1977 

The  second  case  Is  the  storm  of  April  19,1977.  For  this  storm, 
the  southward  neutral  wind  disturbance,  In  figure  23.,  was  already  In 
progress  above  Salnt-Santln  when  observations  started.  Hence,  the 
exact  time  delay  between  the  onset  of  auroral  heating  (vls-a-vls  the 
AE  Index)  and  the  onset  of  the  neutral  wind  disturbance  above  Salnt- 
Santln  cannot  be  determined.  The  observed  neutral  winds  are  generally 
southward  and  show  considerable  structure.  The  model  winds  have  the 
correct  direction,  but  do  not  reproduce  the  observed  structure.  The 
model  clearly  overestimates  the  magnitude  of  the  southward  wind  at  150 
Km. 

For  this  storm,  measurements  of  the  electric  field  are  also 
available.  The  qulet-tlme  electric  field  model  of  Blanc  and  Amsyenc 


<1979)  was  subtracted  from  the  measurements.  Figure  24.  compares  the 
observed  northward  electric  field  due  to  the  disturbance  with  Model 
predictions.  It  is  Important  to  note  that  the  electric  field  observa¬ 
tions  are  the  result  of  two  physical  processes:  the  direct  penetration 
of  the  nagnetospherlc  convection  electric  field  (Vasyllunas,  1970, 
1972),  and  the  Ionospheric  disturbance  dynamo  (Blanc  and  Richmond, 
1980).  The  observations  show  a  northward  electric  field  on  the  order 
of  1-2  mV/m  with  much  structure.  This  is  comparable  to  the  qulet-tlme 
amplitude  measured  by  Blanc  and  Aaayenc  (1979).  The  amplitude  and 
variabllty  is  greatest  from  06UT  until  12UT  when  the  auroral  activity 
is  strong  and  the  penetration  of  the  magnetospherlc  electric  fields  Is 
probably  most  important.  After  12UT  when  auroral  activity  has  de¬ 
creased  the  ionospheric  dynamo  should  have  Its  greatest  relative  im¬ 
portance  on  the  electric  field. 

The  model  predicts  the  correct  direction  of  the  electric  field 
but  lacks  the  structure  shown  In  the  observations  and  overpredicts  the 
magnitude.  This  overpredlctlon  Is  likely  related  to  the  overprediction 
of  wind  strengths  at  150  km,  seen  in  figure  23.  The  generation  of  the 
electric  fields  by  the  disturbance  dynamo  occurs  predominately  below 
ISO  km,  where  daytime  conductivities  are  greatest,  and  the  model  may 
be  significantly  overestimating  these  E  region  winds.  Causes  of  this 
differences  may  be  linked  to  the  Inherent  limitations  of  the  two-dimen¬ 
sional  model  and  to  possible  deviations  of  the  Joule  heating  from  the 
assumed  parameterization  based  on  the  AE  Index. 


Comparison  between  observed  northward  electric  field  and 
predicted  electric  field  over  Saint  Santin  on  April  19, 


VII  CONCLUSIONS  &  FUTURE  WORK 


The  results  of  Richmond  and  Matsushita  (1975)  shoved  that 
thermospheric  circulation  due  to  the  Impulsive  heating  associated  with 
a  geomagnetic  storm  Is  Initially  characterized  by  a  period  of  tran¬ 
sient  gravity  waves.  The  computer  simulations  presented  here  confirm 
that  work,  which  Is  not  surprising  since  I  used  essentially  the  same 
model.  Also,  these  computer  simulations  show  the  development  of  a 
mean  circulation  associated  with  storm-time  heat  Inputs.  This  Is  In 
agreement  with  the  observations  of  Hazaudler  et  al.  (1985)  and  the 
modelling  work  of  Blanc  and  Richmond  (1980).  Also  In  agreement  with 
Blanc  and  Richmond,  the  results  here  show  that  the  storm-time  motions 
of  the  thermosphere  drive  an  electrical  current  system  that  Is  opposite 
In  direction  to,  and  on  the  same  order  of  magnitude  as  the  qulet-tlme 
Sq  current  system.  The  current  generated  by  the  generic  storm  was 
approx.  1/5  the  magnitude  of  the  Sq  current  at  44°  colatitude. 

Both,  the  mean  circulation  and  the  waves,  transport  heat  from  high 
to  low  latitudes.  Heat  transport  by  the  mean  circulation  accounts  for 
2/3  of  the  low  latitude  temperature  Increase.  This  again  confirms  the 
work  of  Richmond  and  Matsushita  (1975).  Richmond  (1979b)  shoved  that 
low  latitude  heating  due  to  heat  transport  via  the  mean  circulation 
dominated  over  heating  due  to  the  viscous  dissipation  of  gravity  waves. 
Not  considered  In  that  work,  however,  was  heating  due  to  heat  transport 
by  the  waves  themselves.  This  mechanism  could  be  important,  especially 
during  the  first  few  hours  of  a  geomagnetic  storm.  It  remains  as  a 
topic  for  future  work. 


Additionally,  work  done  with  C.  Mazaudler  and  A.D.  Richmond 


(Mazaudler  et  al.,  1987)  and  presented  here  (chapter  IV)  showed  that 
the  two-dimensional  model  sucessfully  reproduced  the  large-scale  storm¬ 
time  hvdrodynamlcal  and  electrodynamical  features  of  the  thermosphere 
In  one  case.  In  another  case,  the  model  was  less  sucecessful  In  re¬ 
producing  these  features.  The  details  of  these  features  also  need  to 
be  worked  out . 

Richmond  and  coworkers  are  aware  of  the  limitations  of  the  model 
and  are  working  to  Improve  It.  I  hope  to  participate  In  this  work. 
Possible  Improvements  to  the  model  Include  adding  tidal  effects  to  the 
bottom  of  the  model,  and  using  a  better  parameterization  for  the 
auroral  heating  as  the  current  parameterization  Is  only  accurate  by  a 
factor  of  two. 

The  latitudinal  variations  could  be  better  understood  by  comparing 
the  wind  observations  from  two  incoherent  scatter  radars  at  different 
latitudes  (Millstone  Hill  and  Areclbo)  with  model  outputs.  Longitu¬ 
dinal  variations  that  are  normally  absent  from  2d  simulations  also 
need  to  be  addressed.  This  can  be  done  Ideally  by  using  a  3d  model 
like  the  NCAR  TGCM,  or  by  Incorporating  some  sort  of  longitudinal 
variation  Into  the  current  2d  model. 
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List  of  Symbols: 

radius  of  earth, 
geomagnetic  field, 
radial  component  of  2. 
unit  vector  In  the  direction  of  2. 
specific  heat  of  air  per  unit  mass, 
electrostatic  field, 
magnitude  of  electronic  charge, 
arbitrary  operator, 
acceleration  of  gravity,  9.S  m/s2, 
atmospheric  scale  height, 
variation  of  the  H  component  of  2. 
altitude. 

inclination  of  geomagnetic  field  below  horizontal, 
mass  of  electron  and  Ion  (mean), 
pressure . 

pressure  at  80  Km,  1.038  Pa. 
background  heat  source  per  unit  mass  and  time, 
q  upward  heat  flux. 

R  gas  constant  for  air. 

r  ratio  of  Ion  collision  to  gyrofrequencles. 

T  temperature, 

t  time. 

UJ  hemispheric  Joule  heating  rate. 
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horizontal  air  velocity. 


* 

Z  dimensionless  height  coordinate. 

Y  =  Cp/(Cp-R>. 

€  specific  enthalpy  plus  kinetic  energy  density  per  unit  mass 

of  air. 

0  colatitude. 

icB  molecular  coefficient  of  heat  conduction. 

Kt  turbulent  coefficient  of  heat  conduction. 

molecular  coeffieclent  of  viscosity. 
yt  turbulent  coeffieclent  of  viscosity, 

v  lon-neutral  collision  frequency, 

p  air  density. 

a- 1,02  Pedersen  and  Hall  conductivities. 
ap  height  intergrated  Pedersen  conductivity. 

x  horizontal  viscous  stress. 

♦  east  longitude. 

■P  arbitrary  function  of  time. 

ui  angular  rotation  rate  of  the  earth,  7.29X~^  s“* . 


REFERENCES 


Ahn,  B.,  S.I.  Akasofu,  and  Y.  Kamide,  1983,  The  Joule  heat  production 
rate  and  the  particle  Injection  rate  as  a  function  of  the  geomag¬ 
netic  Indices  AE  and  AL,  J .  Geophys  Res . .  88,  A8,  6275-6287. 

Banks,  P. ,  1977,  Observations  of  Joule  and  particle  heating  In  the 
auroral  zone,  J.  Atmos.  Terr.  Phys.,39.  179-193. 

Banks,  P. ,  J.C.  Feder,  and  J.R.  Doupnlk,  1981,  Chatanlka  radar  obser¬ 
vations  relating  to  the  latitudinal  and  local  time  variation  of 
Joule  heating,  J.  Geophvs.  Res. .  87,  A8,  6869-6878. 

Baumjohann,  W. ,  and  Y.  Kamide,  1984,  Hemispherical  Joule  heating  and 
the  AE  Indices,  J.  Geophvs.  Res..  89,  Al,  383-388. 

Blanc,  M. ,  and  P.  Amayenc,  1979,  Seasonal  variations  of  the  Iono¬ 
spheric  ExB  drifts  above  Saint  Santln  on  quiet  days,  J .  Geophys . 
Res. .  84,  A6,  2691-2704. 

Blanc,  M. ,  A.D.  Richmond,  1980,  The  Ionospheric  disturbance  dynamo, 

J .  Geophys .  Res . ■  85,  1669-1686. 

Brekke,  A.,  1976,  Electric  fields,  Joule  and  particle  heating  In  the 
high  latitude  thermosphere,  J.  Atmos.  Terr.  Phys.,  38,  887-895. 

Brekke,  A.,  and  C.L.  Rino,  1978,  High-resolutlon  altitude  profiles  of 
the  auroral  zone  energy  dissipation  due  to  ionospheric  currents, 

J .  Geonhys .  Res . .  83,  2S17-2524. 

Burnside,  R.G. ,  F-A.  Herrero,  J.W.  Meriwether  Jr.,  and  J.G.  Walker, 
1981,  Optical  observations  of  thermospheric  dynamics  at  Areclbo, 

J .  Geophvs .  Res . .  86,  A7,  5532-5540. 

Burnside,  R.G.,  R.A.  Behnke,  and  J.G.  Walker,  1983,  Meridional  neu¬ 
tral  winds  In  the  thermosphere  at  Areclbo :  Simultaneous  Inco¬ 
herent  and  alrglow  observations,  J .  Geophvs .  Res . .  88,  A4,  3181- 
3189. 

Chapman,  S. ,  1931,  The  absorption  and  dissociative  or  Ionizing  effect 
of  monochromatic  radiation  in  an  atmosphere  on  a  rotating  earth, 
Proc.  Phys.  Soc.  (London).  43,  26-45. 

Chapman,  S. ,  and  J.  Bartels,  1940,  Geomagnetism.  Oxford  Unlv.  Press 
(Clarendon),  London. 


94 


Chlmonas,  G. ,  and  C.O.  Hines,  1970,  Atmospheric  gravity  waves  launch¬ 
ed  by  auroral  currents,  Planet.  Space  Scl..  18,  S65-S82. 

Dickinson,  R.E.,  and  J.E.  Geisler,  1968,  Vertical  motion  field  In  the 
middle  thermosphere  from  satellite  drag  densities,  Hon.  Weather 
Rev. .  96,  606-616. 

Dickinson,  R.E.,  E.C.  Ridley,  and  R.G.  Roble,  1975,  Meridional  circu¬ 
lation  in  the  thermosphere,  I.  Equinox  conditions,  J.  Geophys. 
Res. ■  32,  1737-17S4. 

Duboin,  A.L. ,  and  Y.  Kamlde,  1984,  Latitudinal  variations  of  Joule 
heating  due  to  auroral  electrojets,  J.  Geophvs.  Res.. 89.  24S- 
2S1 . 

Francis,  S.H.,  1975,  Global  propagation  of  atmospheric  gravity  waves: 
A  review,  J.  Atmos.  Terr.  Phvs..  37,  1011-10S4. 

Hanson,  W.B.,  and  H.C.  Carlson,  1977,  The  ionosphere,  The  Upper  Atmo¬ 
sphere  and  Magnetosphere,  National  Academy  of  Sciences,  Washing¬ 
ton  D.C. ,  84-101 . 

Hargreaves,  J.K.,  1979,  The  Upper  Atmosphere  and  Solar-Terrestrial 
Relations,  Van  Nostrand  and  Relnhold  Company,  London. 

Harper,  R.M. ,  1973,  Nighttime  meridional  neutral  winds  near  350  Km 
at  low  to  mid-latitudes,  J.  Atmos.  Terr.  Phvs..  35,  2023-2034. 

Hays,  P.B.,  and  R.G.  Roble,  1971,  Direct  observations  of  thermo¬ 
spheric  winds  during  geomagnetic  storms,  J .  Geophvs .  Res . .  76, 
S316-S321 . 


Hays,  P.B.,  J.W.  Meriwether,  and  R.G.  Roble,  1979,  Nighttime  thermo¬ 
spheric  winds  at  high  latitudes,  J .  Geophvs .  Res . .  84,  1905-1913. 


Hedin,  A.E.,  J.E.  Salah,  J.V.  Evans,  C.A.  Reber,  G.P.  Newton,  N.W. 
Spencer,  D.C.  Kayser,  D.  Alcayade,  P. Bauer,  L.L.  Cogger,  and 
J.P.  McClure,  1977,  A  global  thermospheric  model  based  on  mass 
spectrometer  and  incoherent  scatter  data:  MSIS1,  N2  density  and 
temperature,  J.  Geophvs.  Res.,  82,  2139-2147. 


Hedin,  A.E. ,  1983,  A  revised  thermospheric  model  based  on  mass  spec¬ 
trometer  and  incoherent  scatter  data:  HSIS-83,  J.  Geophvs.  Res.. 
88,  10170-10188. 


Hernandez,  G. ,  and  R.G.  Roble,  1976a,  Direct  measurements  of  nighttime 
thermospheric  winds  and  temperatures  1.  Seasonal  variations  during 
geomagnetic  quiet  periods,  J,  Geophys.  Res. .  81,  2065-2074. 


95 


Hernandez,  G. ,  and  R.G.  Roble,  1976b,  Direct  measurements  of  nighttime 
thermospheric  winds  and  temperatures  2.  Geomagnetic  storms,  J . 
Geophys.  Res..  81,  S173-5181. 

Hernandez,  G. ,  and  R.G.  Roble,  1977,  Direct  measurements  of  nighttime 
thermospheric  winds  and  temperatures  3.  Monthly  variation  during 
solar  minimum,  J.  Geophys.  Res.,  82,  550S-5511. 

Hernadez,  G. ,  and  R.G.  Roble,  1978,  Observations  of  large-scale  thermo 
spheric  waves  during  geomagnetic  storms,  J.  Geophys.  Res..  83,  A12 
SS31-SS38 . 

Hernandez,  G. ,  and  R.G.  Roble,  1984,  The  geomagnetic  quiet  nighttime 
thermospheric  wind  pattern  over  Fritz  Peak  Observatory  during 
solar  cycle  minimum  and  maximum,  J .  Geophys .  Res . .  89,  A1 ,  327- 
337. 

Hines,  C.O. ,  1960,  Internal  atmospheric  gravity  waves  at  ionospheric 
heights,  Can.  J.  Phys..  38,  1441-1481. 

Hines,  C.O. ,  1974,  Propagation  velocities  and  speeds  in  ionospheric 
waves:  A  review,  J,  Atmos.  Terr.  Phy3. .  36,  1179-1204. 

Jacchia,  L.G. ,  196S,  Static  diffusion  models  of  the  upper  atmosphere 
with  empirical  temperature  profiles,  Smithsonian  Contrib.  Astro- 
Phys . .  8. 

Jacchia,  L.G. ,  1971,  Revised  static  models  of  the  thermosphere  and 
exosphere  with  empirical  temperature  profiles,  Spec.  Rep.  332. 
Smithson.  Astrophys.  Observ.,  Cambridge,  Mass. 

Kamlde,  Y. ,  A.D.  Richmond,  and  S.  Matsushita,  1981,  Estimation  of 
ionospheric  electric  fields,  Ionospheric  currents,  and  field 
aligned  currents  from  ground  magnetic  records,  J .  Geophys . Res . , 

86,  801-813. 

Kamlde,  Y.,  and  A.D.  Richmond,  1982,  Ionospheric  conductivity  depend¬ 
ence  of  electric  fields  and  currents  estimated  from  ground  mag¬ 
netic  observations,  J.  Geophys.  Res..  87,  8331-8337. 

Kamlde,  Y.,  and  S.I.  Akasofu,  B.H.  Ahn,  W  Baumjohann,  and  J.  Klsabeth, 
1982,  Total  current  of  the  auroral  electrojet  estimated  from  IMS 
Alaska  meridian  chain  of  magnetic  observatories.  Planet.  Space 
Scl .  .  30,  621-625. 

Killeen,  T.L.,  R.G.  Roble,  R.W.  Smith,  N.W.  Spencer,  J.W.  Meriwether, 
Jr.,  D.  Rees,  G.  Hernandez,  P.B.  Hays,  L.L.  Cogger,  D.P.  Sipler, 
M.A.  Blond! ,  and  C.A.  Tepley,  1986,  Mean  neutral  circulation  in 
the  winter  polar  F  region,  J.  Geophys.  Res..  91,  A2,  1633-1649. 


»V.»TiT*W 


Killeen,  T.L.,  1987,  Energetics  and  dynamics  of  the  earth's  thermo¬ 
sphere,  Reviews  of  Geophysics,  25,  3,  433-454. 

Lyons,  L.R.,  1980,  Generation  of  large-scale  regions  of  auroral  cur¬ 
rents,  electric  potential,  and  precipitation  by  the  divergence  of 
the  convection  electric  field,  J.  Geophys.  Res.,  85,  Al,  17-24. 

Lyons,  L.R.,  T.L.  Killeen,  and  R.L.  Waltersheld,  1985,  The  neutral 
wind  "flywheel"  as  a  source  of  quiet-time  ,  polar-cap  currents, 
Geophys.  Res.  Lett.,  12,  2,  101-104. 

Mazaudier,  C. ,  R.  Bernard,  S.V.  Venkateswaran,  198S,  Saint  Santln 

radar  observations  of  lower  thermospheric  storms,  J .  Geophys .  Res . 
90,  3,  2885-2895,  correction  J.  Geophys.  Res.,  90,  7,  668S-6686. 

Mazaudier,  C. ,  A.D.  Richmond,  and  D.  Brinkman,  1987,  On  thermospheric 
winds  produced  by  auroral  heating  during  magnetic  storms  and 
associated  electric  fields,  Annales  Geophvslcae,  submitted. 

Richmond,  A.D.  and  S.  Matsushita,  197S,  Thermospheric  response  to  a 
magnetic  substorm,  J.  Geophys.  Res..  80,  9,  2839-28S0. 

Richmond,  A.D.,  1979a,  Ionospheric  wind  dynamo  theory:  A  review,  J. 
Geomag.  Geoelectr,,  31,  287-310. 

Richmond,  A.D.,  1979b,  Thermospheric  heating  in  a  magnetic  storm: 
Dynamic  transport  of  energy  from  high  to  low  latitudes,  J.  Geo¬ 
phys  .  Res . .  84,  A9,  5259-5266. 

Richmond,  A.D.,  1983,  Thermospheric  dynamics  and  electrodynamics, 
Solar-Terrestrial  Physics,  B.L.  Carovlllano  and  J.M.  Forbes 
<eds. ) ,  523-607 . 

Rlshbeth,  H. ,  and  O.K.  Garrlott,  1969,  Introduction  to  Ionospheric 
Physics.  Academic  Press,  New  York. 

Rlshbeth,  H. ,  1972,  Thermospheric  wind  and  the  F  region:  A  review, 

J.  Atmos.  Terr,  Phys . .  34,  1-47. 

Roble,  R.G.,  B.A.  Emery,  J.E.  Salah,  and  P.B.  Hays,  1974,  Diurnal 
variation  of  the  neutral  thermospheric  wind  determined  from  In¬ 
coherent  scatter  radar  data,  J.  Geophys.  Res.,  79,  2868-2876. 

Roble,  R.G.,  1977,  The  thermosphere,  The  Upper  Atmosphere  and  Magneto¬ 
sphere  ,  National  Academy  of  Science,  Washington  D.C. ,  57-71. 

Roble,  R.G. ,  A.D.  Richmond,  W.L.  Oliver,  and  R.M.  Harper,  1978,  Iono¬ 
spheric  effects  of  a  gravity  wave  launched  by  the  September  18, 
1974,  sudden  commencement,  J.  Geophys.  Res..  83,  A3,  999-1009. 


Salah,  J.E.,  and  J.M.  Holt,  1974,  Midlatitude  thermospheric  winds  from 


r , 


Incoherent  scatter  radar  and  theory,  Radio  Scl . .  9,  301-313. 


Slpler,  D.P.,  and  M.A.  Blondl ,  1979,  Midlatitude  F  region  neutral 
winds  and  temperatures  during  the  geomagnetic  storm  of  March 
26,  1976,  J .  Geophvs .  Res . .  84,  37-40. 

Slpler,  D.P.,  M.A.  Blondl,  and  R.G.  Roble,  1981,  Midlatitude  F-reglon 
neutral  winds  and  temperatures  during  the  priority  regular  world 
day  14  Aug.  1980,  Planet.  Space  Scl..  29,  12,  1367-1372. 

Slpler,  D.P.,  B.B.  Luokkala,  and  M.A.  Blondl,  1982,  Fabry-Perot  deter¬ 
minations  of  mid-latitude  F-reglon  neutral  winds  and  temperatures 
from  197S  to  1979,  Planet.  Space  Scl..  20,  102S-1032. 

Vasyllunas,  V.M. ,  1970,  Mathematical  models  of  magnetospherlc  convec¬ 
tion,  In  Particles  and  Fields  in  the  Magnetosphere,  edited  by  B.M. 
McCormac,  D.  Reldel,  Hingham,  Hass.,  60-70. 

Vasyllunas,  V.M. ,  1972,  The  Interrelationship  of  the  magnetospherlc 
processes,  In  Earth's  Magnetospherlc  Processes,  edited  by  B.M. 
McCormac,  D,  Reldel,  Hingham,  Hass.,  29-40. 

Vestlne,  E.H.,  1960,  The  upper  atmosphere  and  geomagnetism,  in  Physics 
of  the  Upper  Atmosphere.  J.A.  Ratcllffe,  ed. ,  Chapter  10,  Academic 
Press,  New  York. 

Vickrey,  J.F.,  R.R.  Vondrak,  and  S.J.  Matthews,  1982,  Energy  deposit¬ 
ion  by  precipitating  particles  and  Joule  dissipation  in  the  auro¬ 
ral  ionosphere,  J .  Geophvs .  Res . .  87,  5184-5176. 

Waltersheid,  R.L.,  and  D.J.  Boucher,  1984,  A  simple  model  of  the  tran¬ 
sient  response  of  the  thermosphere  to  Impulsive  forcing,  J,  Atmos. 
Scl. .  41,  1062-1072. 


£HJ> 
fiLtn  el) 

TpllS 

ns-r 


